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ABSTRACT 


The  populations  of  erbium  ions  in  the  tetragonal,  cubic,  and 
trigonal  sites  in  CaF2  crystals  were  measured  for  concentrations 
of  .01%,  .03%,  .1%,  and  .3%  erbium  in  CaF2  crystals.  Experimental 
results  show  that  the  relative  population  of  erbium  ions  in  the 
cubic  (quasi-cubic  or  perhaps  aggregate)  site  increased  as  the 
concentration  of  total  erbium  ions  increased  while  the  relative 
site  population  in  the  tetragonal  site  decreased  with  increasing 
erbium  ion  concentration.  The  relative  population  in  the  tri- 
gonal site  increased  slightly  from  the  .01%  to  the  .1%  concen- 
tration, then  decreased  in  the  .3%  concentration. 

Electron  spin  resonance  techniques  were  used  to  obtain  the 
absorption  resonance  characteristics  of  CaF2:Er^+  in  the 
different  sites.  The  relative  populations  were  determined  as  a 
function  of  total  concentration  of  erbium  ions.  The  doped  crystals 
were  maintained  at  liquid  helium  temperature  (4.2K)  and  spin 
resonance  spectra  taken  as  a function  of  azimuthal  angle  of  the 
DC  magnetic  field  with  respect  to  the  [110]  axis  of  the  crystal  in 
order  to  determine  the  symmetry  and  thus  identify  the  site. 

The  data  was  then  recorded  for  selected  lines,  digitized  and 
processed  in  a PDP-8  minicomputer.  The  data  was  sent  into  a file 
in  the  Honeywell  635  computer  and  received  additional  processing 
using  the  APL  language  in  order  to  arrive  at  the  site  population 
figures. 
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This  project  was  made  possible  by  questions  raised  in  Dr. 

Fontanel  la 's^  research  with  erbium-doped  rare  earth  crystals.  He  has 
found  evidence  that  the  dielectric  constant  for  these  crystals  can  be 

made  temperature  independent  over  a certain  temperature  range  for 

3+  A 

some  particular  concentration  of  Er  ions.  The  real  and  imaginary 

parts  of  the  complex  dielectric  constants  are  given  by  the  Debye 

equations: 

e'  = eH"  + A/[T(1+<A2)] 

and 

e"  = Awt/[T(1+ii>2t2)] 

where  e'  and  e"  are,  respectively,  the  real  and  imaginary  parts  of  the 
complex  dielectric  constant,  e1^  is  the  "high"  frequency  limit  and  A 
represents  the  strength  of  the  dipole.  A = Np2/3e0k  where  N is  the 
concentration  of  a particular  dipole,  p is  the  particular  dipole 
moment,  and  gq  is  the  permittivity  of  free  space.  The  ratio  Np2/T 
can  be  factored  out  of  the  dielectric  constant  equations,  and  for 
some  concentration  of  a particular  dipole  a cancellation  of  the  1/T 
temperature  dependence  can  be  achieved.  Each  of  the  kinds  of  dipoles 
in  erbium-doped  rare  earth  crystals  can  be  associated  with  a parti- 
cular  site  or  local  environment  for  the  erbium  ion  in  the  rare  earth 
crystal  lattice.  Electron  spin  resonance  is  an  independent  technique 
for  distinguishing  the  different  sites  from  one  another  and  also  for 
determining  the  relative  number  of  ions  in  each  site  as  a function  of 
total  concentration  of  erbium  ions  in  the  crystal.  With  this  infor- 
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mation  Dr.  Fontanel  la  believes  that  it  is  possible  to  develop  a tem- 
perature independent  strain  gauge  or  pressure  gauge  which  would  be 
extremely  useful  in  many  military  applications.  In  his  paper1  on 
erbium  doped  CaF^,  Dr.  Fontanella  alluded  to  a site  that  had  not  been 
reported  before.  He  observed  five  different  dipoles  in  his  research, 

but  he  could  find  correspondence  with  only  four  of  the  five  dipoles 

? 

previously  reported.  The  fifth  in  his  work  was  observed  at  low 
temperatures  and  at  high  concentrations  of  erbium  in  the  crystal. 

Since  this  dipole  occurred  at  high  concentrations  and  the  strength 
increased  monotoriical  ly  with  nominal  concentration  he  concluded  that 
this  dipole  was  cluster  associated.  A paper  soon  to  be  published  by 

3 

C.R.A.  Catlow  gives  models  and  arguments  for  an  aggregation  mechanism 
becoming  predominate  as  the  concentration  of  erbium  ions  increases. 
This  aggregate  site  is  said  to  be  associated  with  quasi-cubic  site 

3,4  3+ 

by  many  researchers.’  My  work  with  CaF2:E  shows  that  the  g 
splitting  value  associated  with  this  aggregate  site  is  constant 
for  all  rotational  orientations  about  the  [110]  axis.  This  gives  a 
strong  suggestion  that  the  aggregate  site  geometry  is  at  least 
approximately  cubic. 

By  finding  the  area  of  energy  absorption  peaks  associated  with 
the  different  sites  in  CaF2:Er  f I can  give  Dr.  Fontanella  the 
relative  population  ratios  he  needs  to  explain  the  temperature  com- 
pensation mechanism  in  his  crystals. 
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THEORY  OF  ELECTRON  PARAMAGETIC  RESONANCE 
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The  problem  of  trying  to  determine  the  abundance  of  ions  in  a 
crystal  is  very  basic.  Once  the  ions  are  in  the  crystal,  how  can 
their  presence  be  detected?  For  the  erbium  used  in  my  crystals  the 
detection  technique  is  paramagnetism.  Paramagnetism  is  defined  for 
the  transition  elements  and  is  said  to  exist  whenever  a system  of 
charges  has  a resultant  angular  momentum.1  This  resultant  angular 
momentum  appears  in  the  rare  earth  group  as  ions  having  partially 
filled  inner  electron  shells. 

The  electron  shell  conf iguration  for  the  erbium  atom  is: 


ls22s22p63s23p64s23d104p65s24d105p66s25d04f12. 

In  the  ionic  condition,  two  electrons  are  removed  from  the  6s  shell 
and  one  from  the  4f  shell  giving  a shell  configuration  as  follows: 

1 s22s22p63s23p64s23dT04p65s24d105p66s05d°4f1 1 . 

The  unpaired  electron  in  the  4f  shell  is  shielded  from  the  exterior 
environment  by  the  full  5s  and  5p  shell.  Associated  with  the  angular 
momentum  is  an  intrinsic  magnetic  moment  since  the  electron  can  be 
thought  of  as  a spinning  negative  charge.  (See  Figure  1)*  A free 
electron  has  angular  momentum  L$  = 11/3/4  for  spin  quantum  number 


m_  = 


= ±j-.  An  external  magnetic  field  will  produce  a torque  on  the 

magnetic  moment  causing  the  angular  momentum  vector  to  precess  at 
2 

the  Larmor  frequency  (2ir2.8026  MHz/Gauss  for  a free  electron). 

This  quantizes  the  resultant  angular  moment  L , = mfi,  (See  Figure  2) 
The  magnetic  moment  of  the  unpaired  electron  thus  provides  a handle  by 
which  we  can  grab  onto  the  erbium  ion  and  detect  its  presence.  The 

3 

resultant  magnetic  moment  resolved  along  the  direction  of  the  exter- 
nal magnetic  field  B can  be  written  as: 
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Electron  spin  angular 
momentum 


Assocuated  magnetic  field 


Electron  spin  magnetic 
moment 


For  electrons: 
s = 1/2 

Ls=  /TsTs+TT  = /”3/4 
~^  = 2(1.001159615% 

e = electron  charge 
m * electron  mass 


Figure  1.  Visualization  of  Electron  Angular  Momentum  and 


Magnetic  Moment 


EXTefcNAL 


The  electron  In  an  external  magnetic  field 
has  a resultant  magnetic  moment. 


Figure  2.  Quantization  of  Angular  Momentum 
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u - qiig  ni$ 

where  is  the  Bohr  inannetron 

where  g is  the  spectroscopic  splitting  factor.  For  a free  electron 
q = 2.0023.  This  factor  will  be  very  important  later  on  in  this 
discussion. 

If  an  unpaired  electron  with  energy  E^  is  placed  in  a magnetic 
field  it  experiences  a change  in  energy  AE  = -uB  = gugB.  The 
energy  level  has  split  into  two  levels,  E = En  ± gyDB.  See  Figure  3. 
The  electron  has  the  least  energy  when  the  resultant  magnetic  moment 
is  parallel  to  the  magnetic  field.  If  our  alternating  field  of  fre- 
quency v is  applied  at  right  angles  to  the  external  magnetic  field, 
the  electron  may  undergo  an  energy  level  transition.  This  means  that 
the  electron  absorbed  energy,  hv , and  the  magnetic  moment  flipped 
to  the  anti -paral 1 el  direction.  This  phenomena  is  called  electron 
spin  resonance  or  electron  paramagnetic  resonance.  The  resonance 
condition  is  given  by  aE  = hv  = gugB  and  the  energy  is  made  available 
in  the  form  of  microwaves.  Figure  4 shows  a typical  ESR  set  up. 

The  klystron  is  the  source  of  microwaves  and  the  DC  magnet  pro- 
vides the  external  magnetic  field.  The  cavity  contains  the  sample 
with  the  magnetic  moment.  It  is  one  arm  of  a microwave  bridge  used 
for  detection  of  magnetic  resonance.  Energy  is  absorbed  when  a 
paramagnetic  ion  undergoes  magnetic  resonance.  The  loss  of  energy  is 
detected  as  a leakage  current  in  the  crystal  detector.  This  leakage 
must  be  constant  until  a resonance  occurs  or  false  energy  absorption 
information  will  be  processed  by  the  ESR  equipment. 

The  opposite  arm  from  the  klystron  contains  the  reference  cavity. 
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AE  = -yB  = m guBB  = +gu  B 


Figure  3.  Energy  level  splitting  in  a static  magnetic  field 
showing  the  two  possible  orientations  of  the 


resultant  magnetic  moment. 


Figure  4.  Diagram  of  ESR  System.  Note  the 
klystron,  microwave  bridge(hybrid  "T"),  DC  magnet 
sweep  coils  that  produce  the  first  derivative  output 
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Figure  4a.  DC  magnet  and  K-band  klystron  in 

room  temperature  configuration.  Note : a 1 1 experiments 

in  this  report  used  X-band  at  liquid  helium  tempera- 
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This  is  necessary  in  order  to  ensure  standing  waves  in  the  waveguide 
and  to  balance  the  bridge.  (See  Appendix  B for  details  on  using 
the  ESR  equipment.)  The  ESR  instrumentation  takes  the  first  derivative 
of  the  absorption  curve.  This  locates  the  resonance  point  much  more 
accurately.  However  if  the  Phase  Tuning  on  the  klystron  is  not 
correct  the  first  derivative  output  will  be  asymmetrical  with 
respect  to  the  base  line  of  the  signal.  (See  Appendix  E*fo r discussion 
of  distortion.)  This  distortion  shows  up  later  when  the  absorption 
curve  is  recovered.  If  a free  ion  with  an  unpaired  electron  is 
placed  in  a magnetic  field,  the  ion  energy  levels  are  given  by 
E = Eg  ± gvigB  where 

g = 1 + J(J+1)+S(S-H)-L(L+1)  ^ 

2J( J+l ) 

and  L,  S,  and  J,  are  the  orbital,  spin,  and  resultant  angular 
momenta,  respectively.  The  extra  frills  on  g are  to  take  into  account 
all  the  other  electrons  in  the  ion.  The  ion  can  be  said  to  have  a 
resultant  magnetic  moment. 

If,  now,  the  erbium  ion  is  placed  in  a crystal,  the  g-factor 

becomes  a little  more  environment-dependent.  The  crystals  used  in 

this  study  are  CaF2*  Calcium  is  a divalent  cation,  while  erbium  is 

a trivalent  cation.  Suppose  a calcium  ion  is  replaced  with  an  erbium 

ion.  There  must  be  a charge  compensating  ion  somewhere  in  the  crystal 

to  keep  the  crystal  electrically  neutral.  The  dominant  charge  com- 

3+ 

pensating  mechanism  for  CaF2:Er  is  an  interstitial  fluorine  ion. 

The  crystals  were  grown  in  a fluorine  atmosphere  and  are  very  pure. 

2 4 

Thus  the  (j  and  0H~  compensation  mentioned  in  Heist  and  Fongq  (1970) 

*see  apl  Pftot<xfyv\  ••()•*» Tuftt" 
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is  not  a factor.  The  position  of  the  charge  compensating  fluorine 
ion,  with  respect  to  the  erbium  ion,  will  affect  the  local  geometry 
of  the  erbium  ion.  This  means  that  the  g-splitting  factor  is  now 
a tensor,  and  is  dependent  on  spatial  orientation, 

In  CaF^  the  fluorine  ions  form  a cubic  lattice  structure  and  the 
calcium  ions  form  a face-center  cubic.  Each  calcium  is  surrounded 
by  eight  fluorine  ions  and  are  in  alternating  fluorine  cube  centers. 

In  other  words  alternating  fluorine  cube  centers  are  empty  and  can 
be  filled  with  a compensating  fluorine  ion.  If  the  fluorine  occupies 
the  center  of  fluorine  ions  in  the  nearest  neighbor  site  along  a 
[100]  axis  from  the  erbium  ion,  the  erbium  ion  sees  4-fold  or 
tetragonal  symmetry.  (See  Figure  5).  If  the  charge  compensating 
fluorine  is  in  the  next-nearest  neighbor  site  along  a [111]  axis 
from  the  erbium  ion,  the  erbium  ions  see  3-fold  or  trigonal  symmetry. 

If  the  fluorine  is  very  far  from  the  erbium  ion  the  local  symmetry 
of  the  erbium  ion  is  not  distorted  and  the  erbium  ion  sees  cubic 
symmetry. 

Figure  6 shows  the  possible  fluorine  locations  for  tetragonal 
and  trigonal  sites  although  there  are  six  tetragonal  sites,  they  are 
degenerate  in  pairs  along  [001],  [010],  and  [100]  axes,  so  there 
appears  to  be  only  3 equivalent  sites.  There  are  eight  trigonal  sites, 
but  they  are  degenerate  in  pairs  along  [111],  [HIT,  [111],  and  [Til] 
directions,  so  there  appears  to  be  only  four  equivalent  trigonal  sites. 
The  tetragonal  sites  can  be  distinguished  from  the  trigonal  sites 
by  choosing  the  proper  rotational  axis  of  the  crystal.  As  mentioned 
earlier,  the  g-factor  is  now  a tensor  and  depends  on  the  spatial 
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orientation  of  the  crystal  with  respect  to  the  external  magnetic 
field.  In  Figure  7 the  general  case  of  g as  a function  of  rotation 
angle  is  shown.  For  some  rotational  angles,  g defines  the  change  of 
the  energy  level  splitting  as  the  external  magnetic  field  increases, 
that  is 


As  the  crystal  is  rotated  in  the  magnetic  field  about  a given  axis, 
g changes  as  illustrated  in  Figure  7.  The  set  of  all 

g's  for  all  rotations  about  all  axes  is  defined  as  the  splitting 
ellipsoid.  The  distance  from  the  center  of  the  ellipsoid  (the 
erbium  ion)  to  the  surface  of  the  ellipsoid  in  the  direction  of  the 
external  magnetic  field  gives  the  value  of  g.  Note  that  as  g changes 
the  AE  = hv  from  the  microwaves  remains  constant.  The  frequency  of 
the  microwaves,  v,  is  constant  therefore  giving  constant  energy  hv  for 
any  resonance  that  take  place.  The  magnetic  field  (in  Gauss)  is 
varied  over  a certain  range  and  whenever  the  resonance  condition 
is  satisfied,  aE  = hv  = 9PgBext>  a resonance  occurs.  As  the  crystal 
is  rotated, g changes  and  the  resonance  occurs  at  a different 
magnetic  field  strength  for  that  orientation.  Figure  8 shows  the 
change  in  g as  the  ellipsoid  is  rotated  in  the  magnetic  field. 

5 

C.  W.  Rector  analyzed  several  rotational  axis  and  found  that  the 
[110]  axis  is  best  for  distinguishing  the  tetragonal  sites  from  the 
trigonal s.  In  the  rotational  spectrum  of  g values,  two  tetragonal 
sites  show  up  as  one  line  and  the  third  appear  in  a second  line. 

Three  trigonal  lines  appear  in  the  spectrum,  two  behaving  the  same 
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Figure  7.  Relationship  of  Energy  Level  Splitting  and 
g-factor.  As  the  crystal  is  rotated  about  a synmetrv 
axis  in  an  external  magnetic  field  the  g-splitting  factor 
changes.  Note:  Since  v is  constant,  the  resonance  energy 

is  constant  but  the  resonant  magnetic  field  strenth,  Bp  , 
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way,  thus  appearing  as  one  line,  giving  only  3 lines  instead  of  four. 
Figure  9 shows  the  g spectrum  for  the  sites  as  a function  of  rotational 
angle.  The  equations  for  the  g of  each  line  are  given  below.  A 
program  "ESRPLOT"  is  given  in  the  Appendix  E to  reproduce  the  g-map. 

A listing  of  the  values  at  5°  intervals  may  also  be  found. 
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METHOD  OF  RESEARCH 


The  first  step  in  gathering  resonance  absorption  data  is  to  orient 
the  crystal  in  the  microwave  cavity  with  the  [110]  axis  perpendicular 
to  the  magnetic  field.  The  procedure  for  this  is  given  in  Appendix  A. 
The  [110]  axis  becomes  the  rotational  axis  for  the  external  DC  magnet. 
Once  the  crystal  is  properly  mounted,  the  magnetic  resonance  road  map 
(See  Figure  10)  is  used  to  locate  0°  or  the  [001]  azimuthal  orienta- 
tion. I found  it  more  useful  to  map  the  magnetic  resonance  points 
as  a function  of  rotational  angle  rather  than  the  g splitting  values. 
The  map  is  made  from  the  g splitting  equation  by  this  equation: 


B - -h-v~. 


where  v is  the  frequency  of  the  microwaves.  The  frequency  v is  deter- 
mined by  using  a small  amount  (a  few  grains)  of  DPPH  taped  to  the 
side  of  the  cavity  DPPH  has  a fixed  g value  of  2.0036  ± .0003  and  v 


can  be  found  by  this  formula: 


(2.0036)„b  BMt 


where  is  the  magnetic  field  strength  for  which  DPPH  will 

resonant.  (C  t = 3.29  to  3.5  K9)  Figure  10  gives  a graph  of  the 
B spectrum. 

I found  that  the  tetragonal  and  cubic  resonance  peaks  have  a 
greater  separation  at  0°  (or  180°).  Other  orientations  that  looked 
promising  on  the  road  map  proved  to  have  too  many  hyperfine  or 
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Figure  10.  Magnetic  Road  Map 
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harmonic  resonances  that  crowded  the  resonance  spectrum.  ideally  I 
wanted  to  find  an  isolated  peak  corresponding  to  a particular  site. 
However,  1 found  that  the  peaks  overlap.  Hie  best  I could  do  was  to 
minimize  peak  overlap. 

The  data  was  taken  at  liquid  helium  temperatures  (4.2°K)  to 
reduce  lattice  vibrations  and  to  reduce  the  spin-orbit  interaction, 
which  broadens  the  absorption  curves,  making  an  accurate  separation 
of  the  individual  curves  very  difficult.  Working  at  4.2°K  the 
lattice  vibrations  are  "frozen  out"  and  the  absorption  peaks  become 
sharper  and  more  distinct,  producing  better  peak  separation. 

The  trigonal  site  can  be  isolated  since  its  resonance  spectrum 
has  such  a wide  variation  in  magnetic  field  strength.  I used  the 
resonance  ( 1.74KG)  at  an  azimuthal  rotation  of  +39°  from  the  angle 
which  I used  for  the  tetragonal  sites.  Since  the  microwave  cavity 
is  not  cylindrical  it  presents  a different  cross-section  to  the 
magnetic  field  and  as  the  the  magnet  is  rotated  from  0°  magnet  to 
60°  magnet,  the  signal  amplitude  decreases.  DPPH  is  used  as  a 
reference  absorption  curve  for  the  two  orientations  since  its  g value 
does  not  change  with  rotation  in  a magnetic  field. 

The  identification  of  the  sites  at  0°  (hereafter  called  9°-theory 
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and, when  in  reference  to  magnet  orientation,e°-  mag)  is  shown  in 
Figure  11.  The  larger  tetragonal  peak  accounts  for  2/3  of  the  total 
tetragonal  sites,  while  the  smaller  peak  accounts  for  the  other  third. 
For  the  trigonal  spectrum  the  two  lines  that  vary  from  g = 8.55  to 
g = 3.29  each  represent  1/4  of  the  total  trigonal  site  population. 

The  third  trigonal  line  accounts  for  2/4  of  the  total  trigonal  site 
population. 

With  the  sites  identified  I began  the  real  work  of  getting  areas 
under  absorption  curves.  The  basic  problem  I encountered  here  was 
getting  a consistent  symmetric  line  shape  in  the  first  derivative 
curve.  Appendix  B gives  the  details  for  handling  the  equipment  for 
best  line  shape.  Appendix  Z has  a discussion  about  the  software 
method  I used  to  correct  for  phase  distortion  in  the  first  derivative 
output.  I also  proposed  a few  variations  on  the  method. 

After  the  resonance  signal  is  as  symmetric  as  possible  the  signal  is 
recorded  on  a Hewlett  Packard  instrumentation  tape  recorder. 

This  is  also  covered  in  Appendix  B.  As  soon  as  all  the  recordings  are 
done,  they  are  taken  to  a PDP-8E  minicomputer  which  has  analog  to 
digital  capabilities.  The  analog  voltage  signal  recorded  on  the  tape 
is  converted  to  numbers  stored  in  the  PDP-3  memory.  The  procedure 
for  this  and  for  sending  the  numbers  to  the  DTSS  (Dartsmouth  Time 
Sharing  System)  system  are  in  Appendix  C.  One  item  of  interest  is  the 
low  pass  moving  average  digital  filter  incorporated  in  the  software. 
This  is  discussed  in  Appendix  G.  The  result  of  this  filter  gives 
outstanding  "smoothness"  to  the  data.  The  program  is  constructed  so 
that  either  the  filtered  or  non-filtered  data  may  be  displayed  on  an 


oscilloscope.  Figure  12  gives  an  example  of  a DPPH  curve  with  and 
without  filtering.  The  operator  has  the  option  of  controlling  the 
filter  parameter  which  is  the  width  of  the  moving  summation  window. 

Once  the  data  is  transferred  to  DTSS,  several  APL  programs  are  used 
to  manipulate  the  data  to  recover  the  absorption  curve  and  get  the 
area  under  it. 

By  following  Appendices  A thru  E a researcher  can  duplicate  the 
procedure  I used  for  the  experiment.  The  following  is  a summary  of 
the  action  of  each  APL  program, 

RAWCOR:  This  program  eliminates  the  PDP-8  modulo  409G  wrap-around, 

removes  DC  bias,  and  corrects  for  initial  potato.  Result 
is  called  COROAT. 

CUTPLOT:  Is  used  to  cut  down  the  number  of  points  to  be  plotted. 

In  this  way  a file  of  4096  data  points  may  be  reduced  to 
200  plotting  points.  See  Appendix  D for  use  of  plotted 
graph.  Result  is  called  CPOUT. 

DISTORT:  This  program  corrects  for  phase  distortion  produced  by 

inbalance  in  phase  tuning  on  klystron.  It  assumes  a 
Gaussian  error  in  the  leakage  current  that  detects  the 
magnetic  resonance.  Result  is  called  CORD. 

INTTEG:  This  program  performs  indefinite  integration  on  the  first 

derivative  data.  The  range  of  integration  is  an  input  to 
the  program.  The  result  is  called  ICURV.  The  output  is 
the  absorption  curve  with  an  integration  ramp. 
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ABSCOR: 
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FIT:  This  program  fits  a particular  absorption  curve  to  a 

Gaussian  curve,  The  boundaries  of  the  curve  are  inputs. 
Amplitude  A and  exponential  coefficient  B are  outputs;  ST 
is  standard  deviation.  The  mid-point  of  the  Gaussian 
curve  is  found  by  the  program.  The  result  is  called  CURFIT. 
INTSIM:  This  program  performs  the  definite  integration  on  the 

fitted  curve  to  arrive  at  a number  equal  to  the  area  under 
the  fitted  curve. 

APLTOBAS:  This  program  is  necessary  for  converting  APL  data  files 

into  standard  graphic  format.  Use  CUTPLOT  while  in  APL; 

exit  from  APL;  run  APLTOBAS  (a  BASIC  program). 

ESR  PLOT:  This  program  outputs  the  g-splitting  or  magnetic 

3+ 

resonance  point  spectra  for  CaF,,:Er  in  table  form  or 
in  a form  suitable  for  plotting,  (a  BASIC  program).  It 
requires  "Split"  as  an  output  file,  (from  -30  to  390°) 

Figure  14  shows  a pictoral  procedure  for  recovering  the  absorp- 
tion curves.  Figure  14a.  shows  the  phase  distortion  of  a first 
derivative  curve  (dark  line)  and  its  corrected  shape.  This  is  done 
with  the  program  called  "DISTORT."  Figure  14b.  shows  the  DPPH  first 
derivative  line  and  14c.  shows  the  indefinite  integration  results  (INTTEG). 
Figure  15  shows  the  results  of  a complex  erbium  spectrum  and  shows 
how  the  program  "FIT"  will  separate  the  curves.  All  that  remains  is 
to  use  "INTSIM"  to  find  the  area  of  the  fitted  curves. 


Figure  12.  Unfiltered 
DPPH  resonance 
signal . 


Figure  13.  results  of  digital  filter  on  17  Hz 
beat  frequency  in  ESPv  equipment. 
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RESULTS 

I was  able  to  identify  the  tetragonal,  cubic,  and  trigonal 
sites  in  the  magnetic  resonance  spectrum  of  erbium  doped  CaF2  using 
the  magnetic  resonance  raad  map.  Figure  16  shows  an  absorption 
curve  for  .1%  erbium  concentration  with  its  associated  sites.  This 
was  taken  at  0°  azimuthal  angle  with  respect  to  the  [110]  axis  of 
rotation.  The  values  of  the  magnetic  field  strength  for  these  peaks 
are  given  in  the  Appendix  G. 

I was  not  always  able  to  identify  the  trigonal  site 
between  the  tetragonal  and  the  cubic  sites  as  shown  in  Figure  16. 

I would  then  rotate  to  +39°  to  observe  the  trigonal  site  at  about 
1.74  kilogauss  and  record  it  separately  using  the  DPPH  resonance 
as  the  reference. 

In  the  .3%  concentration  I found  the  cubic  site  becoming 
dominant  in  absorption.  I noted  that  in  the  rotational  study  of 
that  concentration,  the  g-splittina  value  associated  with  the 
cubic  did  not  change.  As  I stated  in  the  introduction  to  this 
paper,  this  cubic  site  has  recently  been  thought  to  be  an  aggregate 
or  cluster  of  erbium  ions  and  the  associated  interstitial  fluorines, 
fty  studies  show  that  the  cluster  must  be  at  least  cubic  in  nature. 

My  results  on  the  relative  populations  of  the  sites  as  a 
function  of  total  concentration  of  erbium  ions  is  shown  in  Figure  17. 

The  cubic  becomes  dominant  as  the  concentration  of  erbium  increases. 
Perhaps  the  tetragonals  are  crowding  eachother  and  becoming  aggregates. 
Figures  18  through  23  shows  some  first  derivative  curves  and  t.heir 
associated  absorption  curves.  In  Figure  23,  the  .OH  erbium,  I 
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had  a little  extra  to^dohn  order  to  recover  the  absorption  curve 
that  would  be  suitable  for  definite  integration.  I took  the  data 
for  the  .01%  before  I realized  the  serious  effects  of  phase 
distortion.  Figure  24  gives  the  relative  percent  figures  for  each 
concentration  of  erbium. 
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Figure  16. 
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APPENDIX  A 


PROCEDURE  FOR  MOUNTING  A F.C.C.  CRYSTAL 

1.  Orient  GYMBAL  JIG  for  perpendicularity.  Make  sure  the  crystal  end 
of  the  mounting  rod  maintains  a fixed  location  in  space  for  any  rota- 
tion . 

2.  Orient  laser  to  shine  across  diameter  of  Spencer  Spectrometer 
Table  at  a height  just  above  mounting  platform  (Silica  Glass)  lying  on 
the  table  (Fioure  A1 ) . 

3.  Make  sure  table  is  level  by  using  a prism  to  reflect  the  beam 

360"  about  table  onto  chart  paper.  If  bean  hits  paper  at  equal  distance 
from  bottom  of  paper,  through  360°;  then  table  is  level. 

4.  Wash  all  gear  with  acetone  and  rinse  with  methyl  alcohol. 

5.  For  F.C.C.  crystals  cleavage  planes  are  [111]  planes  two  [111] 
planes  intersect  in  [110]  line.  Therefore  you  must  get  two  good  cleav- 
age planes  to  intersect  forming  a [110]  axis  (Figure  A2). 

C.  Use  full  strength  duco  cement  and  glue  the  mounting  rod  to  the 
crystal  as  close  to  the  [110]  axis  as  you  can  (Figure  A3). 

7.  Insert  rod  into  teflon  swivel  (white  and  adjust  length  to  be 
just  above  the  silica  glass.  Put  rubber  cement  on  outer  surface  of 
swivel  (Fiqure  A4) . 

8.  Orient  crystal  so  that  laser  beam  is  reflected  to  same  height 
above  table  for  both  faces.  To  compare  one  face  to  another,  move 
laser  thru  small  angle.  Rotate  table  and  thus  crystal  to  see  if 
crystal  reflection  height  is  constant  (Figure  45). 

9.  Prepare  dilute  duco  with  40,v  acetone  added. 


A2 


10.  Align  silica  qlass  as  shown  in  Figure  A7. 

11.  Apply  dilute  duco  to  NON  used  surfaces.  Continue  to  monitor 
reflections  until  glue  is  dry.  The  reflections  should  not  move  for 
proper  alignment  of  [110]  axis  arctan(Y/X)  < 3°  (Figure  A6). 

12.  When  dry  (an  hour  at  least)  wrap  a layer  of  rubber  cement  around 
crystal  and  top  of  silica  glass.  This  holds  crystal  in  place  at 
liquid  helium  temperatures  in  case  the  duco  becomes  brittle  and  breaks 
its  bond. 

13.  When  dry  (30  minutes)  orient  crystal  in  bottom  of  microwave  cavity 
as  in  Figure  A7.  Thus,  0°  location  on  road  map  will  nearly  be  0° 
magnetic  when  you  start  looking  for  axial  orientation  of  crystal. 

See  road  map  Appendix  B.  Put  small  amounts  of  rubber  cement  on  sides 
of  silica  glass  but  not  under  it  (Figure  A3).  Allow  30  minutes  to 
dry. 

14.  For  reference  & sample  of  DPPH  should  be  placed  in  the 

cavity.  Open  the  DPPH  bottle  and  insert  an  extended  paper  clip  all  the 
way  to  the  bottom  of  the  bottle.  Some  DPPH  will  adhere  to  the  wire, 
and  this  is  a sufficient  amount  of  DPPH  to  use.  Rub  off  some  of  the 
DPPH  onto  a piece  of  tape  and  cut  tape  away  from  rest  of  tape.  Attach 
tape  to  wall  of  cavity. 

15.  Attach  cavity  to  wave  guide  as  in  Figure  A9.  You  are  now  ready  to 
go  to  Appendix  B. 


Figure  A5. 
Laser. 
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Figure  A6.  Mounting  Crystal  to  FAC6 
Silica  glass. 
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Figure  A8.  Attaching 
Crystal  to  Cavity  Floor 
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Figure  A7.  Orientation  of  Crystal  in  Cavity 
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Figure  A9.  Attaching  Cavity  to  Waveguide 
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APPENDIX  B 


B1 


SET-UP  OF  ESR  EQUIPMENT 

1.  Turn  on  cooling  water  - next  to  NMR  magnet  power  supply. 

NOTE:  90psi  before  filter. 

80  psi  after  filter 
20  psi  on  magnet 

water  flow  into  large  clyindrical  drain 

water  flow  into  small  cylindrical  drain  for  your  particular 
KLYSTRON  (X  or  KBAND) 

(For  my  work  I used  X band.) 

2.  Turn  on  power  to  Westinghouse  "Main  A"  box  in  corner  on  left  side 
of  magnet  bench  (up  is  on). 

3.  Make  sure  fan  in  'Feldial  Mark  II"  panel  in  varian  consol  is  rotat- 
ing. You  can  hear  it  vibrating  if  it  is  working.  Failure  of  this  fan 
causes  "Potato'^"  An  uncontrollable  drift  of  voltage  output  to  strip 
recorder  pen  to  the  left  or  right.  (The  fan  cools  power  transistors 

in  the  circuitry. ) 

4.  Follow  starting  instructions  in  varian  manual  V-4502  Operating 
Instructions  PUB  number  87-114-200  D368. 

High  Power:  Page  3-5  to  3-27 

Low  Power:  Page  3-29  to  3-56 

Since  I used  low  power  with  a reference  cavity  all  my  comments  will 
apply  to  low  power  only. 

a.  Start  on  page  3-31  read  precaustions . 

b.  Initial  setup  ends  on  page  3-34.  Let  system  warm  up  at  least 
3 hours  and  up  to  12  hours. 

5.  While  systewis  warming  up  refer  to  Procedures  for  Handling  Liquid 
He! ium  PUB  Number  87-250-031  A865.  Page  14-19  "Pre  Cooling  Procedure" 
and  page  1-1  to  1-7  for  various  precautions  and  interesting  facts 


about  liquid  helium.  It  is  recommended  that  gloves  and  safety  mask 
be  used  while  working  with  liquid  nitrogen  and  helium. 

6.  Open  wooden  cabinet  containing  waveguide  and  cryostats. 

7.  Carefully  lower  waveguide  and  cavity  (See  Appendix  A)  into 
hel ium  cyrostat. 

8.  Carefully  lower  helium  cryostat  into  nitrogen  cryostat.  Be  sure 
that  styrofoam  insulator  is  in  place  around  liquid  helium  cryostat. 

9.  Make  sure  that  the  helium  cryostat  is  not  touching  the  nitrogen 
cryostat  by  observing  that  the  helium  cryostat  is  hanging  free  in  the 
nitrogen  cryostat. 

10.  Slide  the  cryostat  and  waveguide  assembly  on  to  rails  above 
magnet  and  center  the  cavity  between  the  pole  faces  and  helmholtz  coils. 

11.  About  every  10-12  runs  re-establish  the  vacuum  in  the  helium 
cryostat.  Then  leak  a small  amount  of  nitrogpn  gas  into  the  vacuum. 

See  page  14-15  in  Procedures  for  Handling  Liquid  Helium  . 

12.  Connect  the  flexible  waveguide  and  H-bend  section  between  klystron 
and  cryostat.  Be  sure  screws  are  tight. 

13.  Before  starting  liquid  nitrogen  transfer  start  a 5 lb/sec  flow  of 
helium  gas  into  helium  valve  on  top  of  the  microwaveguide.  This 
flushes  out  atmospheric  air  and  moisture  which  would  condense  and  form 
ice.  (This  would  draw  microwave  power  and  give  false  resonance  in- 
formation.) Be  sure  to  open  the  other  helium  valve  and  put  a rubber 
stopper  in  the  remaining  hole  on  top  of  the  waveguide  section. 

14.  Transfer  liquid  nitrogen  to  the  top  of  the  nitrogen  cryostat. 

15.  Go  back  to  V-4502  Op.  Instructions  page  3-35  and  continue  thru 
page  3-37.  This  sets  up  the  dip  in  the  power  mode  of  the  klystron  on 
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the  scope.  (Use  signal  level  = 1600) 

16.  Liquid  helium  transfer: 

a.  Place  pressurizing  tube  on  liquid  helium  dewar.  (See  Figure  B1 ) 

b.  Take  out  rubber  stopper  on  top  of  waveguide  section. 

c.  Put  on  gloves. 

d.  Take  transfer  tube  and  insert  it  into  cryostat  and  liquid 
helium  dewar.  (See  Figure  B2) 

e.  Stop  helium  gas  flow  and  connect  helium  gas  to  pressurizing 

tube. 

f.  Start  5 Ib/sec  flow  of  helium  gas. 

g.  Watch  scope  for  movement  of  dip  as  cold  He  gas  enters  the 
cavity.  It  will  shrink,  increasing  the  frequency,  and  dip  will  move  off 
to  the  right. 

h.  As  liquid  helium  enters  the  cavity  the  dip  will  move  very 
quickly  to  the  left.  This  is  because  of  the  change  of  dielectric 
constant  from  gas  to  liquid  in  the  cavity. 

i.  When  liquid  helium  level  is  one  inch  below  the  level  of  the 
liquid  nitrogen  stop  helium  gas  flow. 

j.  Disconnect  rubber  tube  between  helium  dewar  and  helium  gas 
tank.  This  allows  excess  helium  gas  pressure  to  leak  off. 

k.  Put  on  gloves. 

l.  Remove  transfer  tube  and  replace  it  on  the  wall  mount. 

m.  Put  in  rubber  stopper  and  close  helium  valves  on  top  of 
waveguide. 

17.  Start  again  in  V-4502  Op.  Instructions  page  3-37  (8)  and  continue 
thru  page  3-41.  Be  sure  to  go  to  page  2 after  you  finish  page  3-41. 

This  talks  about  adjusting  the  crystal  detector. 
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18.  Set  up  for  recording  resonance  signal  on  strip  recorder.  Continue 

with  page  3-42  to  3-50.  I used  generally:  Sig  Lev  630 

Swp  Frq  400 
Ref  Pha  3.8  (RED) 

Leakage  250uA 
Response  1 

Sig  Lev  from  2 to  630 

19.  To  adjust  reference  phase: 

a.  Find  OPPH  resonance  (3.34-3.30  KG). 

b.  Set  range  to  250. 

c.  Move  percent  of  range  dial  until  recorder  is  driven  to  a maxi- 
mum swing.  Adjust  referencephase  until  maximum  is  achieved. 

20.  You  are  now  ready  to  find  axial  orientation  of  crystal.  To 
familiarize  yourself  with  the  magnetic  resonance  look  at  Figure  B3. 

This  is  the  road  map.  Foranother  copy  of  the  road  map,  See  Appendix  E 
under  ESRPL0T.  InAppendix  A I told  you  how  to  get  the  0°  or  180° 
orientation  near  the  0°  magnet  orientation.  The  procedure  for  getting 
the  0°  theoretical  is  as  follows: 

a . Set  up  Fieldial  panel  : 

Time  1 min 

Dial  (Field  set-KG)  1000 
Range  1 . 0 KG 
Triangular  Sweep 

b.  Start  at  0°  magnet  and  take  a one  minute  sweep  through  1 KG. 

c.  Rotate  5°  and  repeat  until  you  have  gone  thru  '60°.  This  will 
give  you  an  excellent  view  of  the  resonance  spectra. 

d.  Look  for  similarities  in  resonance  patterns  on  each  side  of 
0°  magnet.  Note  on  road  map  that  spectrum  is  identical  on  either  side 
of  0°  theoretical  and  a crossing  occurs  at  approximately  * 55° . From 
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the  symmetry  of  resonance  spectra  you  should  be  able  to  get  within 
2.5°  of  true  0°  theoretical.  The  whole  idea  is  to  find  0°  theoretical 
because  this  angle  has  the  best  separation  of  the  different  site 
absorption  curves. 

21.  You  are  now  ready  to  worry  about  symmetry  of  the  first  derivative 
curve. 

a.  Read  page  5-11  to  5-14  in  EPR  Operational  Techniques  PUB 
Number  87-114-402  C864  and  page  2 to  3 in  V4502  Operating  Instructions. 

b.  Using  the  following  settings  you  can  eliminate  many  of  the 
problems : 

(1)  On  the  klystron:  attenuation  4.0 

cryostat  isolator  2.0 

(2)  On  the  varian  console:  time  10  minutes 

response  1.0 
sweep  field  630 

Reference  phase  should  not  be  changed  from  step  18. c.  These  settings 
will  correct  saturation,  too  fast  scan  rate,  and  give  correct  modula- 
tion for  line  shape. 

22.  You  are  now  ready  to  try  to  make  a particular  first  derivative 
curve  symmetric. 

23.  a.  For  C F,  doped  with  erbium  the  main  cubic  and  tetragonal 

a l 

resonances  occur  within  500  to  1500  gauss.  Set  magnet  for  0°  theoreti- 
cal . 

b.  Set  TIME  at  10  minutes. 

c.  Set  RANGE  at  1 KG. 


d.  Set  FIELD  SET  dials  at  1.000  KG. 
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e.  Set  POR  DIAL  (PERCENT  OF  RANGE  DIAL)  at  -45  this  should  be  a 
region  of  magnetic  field  in  which  there  are  no  major  resonances. 

f.  Center  the  strip  recorder  pin  on  the  chart  paper. 

g.  Move  the  POR  DIAL  to  +45.  If  the  pen  does  not  return  to  center 
you  have  a major  "Potato"  problem. 

h.  Monitor  AFC  "Cat's  Eye"  on  varian  console  scope.  Probably 
the  AFC  will  be  shifted  to  the  extreme  left  or  right. 

i 

i.  If  pen  shifted  to  the  left  you  have  up  potato;  to  the  right 
is  down  potato.  Down  potato  has  an  AFC  shifted  to  the  left.  Up 
potato  has  an  AFC  shifted  to  the  right. 

j.  Move  the  PHASE  TUNING  on  the  klystron  either  left  or  right 

f ! 

to  move  the  AFC  back  towards  a symmetrical  display. 

k.  To  increase  PHASE  TUNING  you  must  rotate  the  dial  to  a setting 
much  greater  than  the  value  you  want,  then  rotate  it  back  to  the 
desired  value.  Same  for  decreasing  it  an  extreme  rotation  past  the 
desired  value.  Then  a return  to  the  desired  value.  This  is  done 

to  momentarily  confuse  the  AFC  and  make  a change  possible. 

l . After  any  change  to  the  PHASE  TUNING  the  AFC  ERROR  must  be 

corrected  and  LEAKAGE  CURRENT  (Controlled  by  the  SLIDE-SCREW  TUNER 
MICROMETER  on  right  side  of  klystron)  must  be  maintained  at  250uA 
when  not  in  resonance.  NOTE:  Leakage  current  increases  with  a clock- 

wise rotation  of  the  slide-screw  tuner  and  decreases  with  counter- 
clockwise rotation. 

m.  Go  to  e.  until  pen  shift  is  gone. 

n.  Find  a major  peak  in  the  resonance  spectra  and  set  the  POR 
DIAL  for  5 % less  than  that  peak's  resonance  value.  Using  the  TIME  = 

Ll  ' 
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10  MIN  and  RANGE  = lKG^ punch  INCREASE  and  sweep  through  to  the 
resonance.  Now  if  one  half  of  first  derivative  curve  is  shorter  than 
the  other,  you  must  correct  the  phase  tuning  according  to  the  follow- 
ing rules: 

To  shift  symmetry  right  --  increase  phase  tuning. 

To  shift  symmetry  left  --  decrease  phase  tuning. 

The  method  of  change  is  outlined  in  k, 

o.  Go  to  n.  until  you  are  satisfied  with  symmetry.  NOTE:  You 

can  never  get  all  of  the  phase  distortion  out  but  the  more  symmetric 
the  better.  The  program  "distort"  will  take  out  the  rest  of  the  dis- 
tortion, 

p.  (1)  Now  set  up  tape  recorder  as  shown  in  Figure  B4. 

(2)  Recorder  switch  on  EPR  consol  must  be  turned  to  PDP-3 
position. 

(3)  Signal  level  and  the  external  amplifier  must  be  adjusted 
so  as  not  to  saturate  the  external  amplifier  (741  chip 
±9  volts  output  range). 

NOTE:  External  amplifier  has  3 positions: 

X25 

X50 

X500. 


The  signal  has  to  go  thru  one  of  these.  There  is  a switch  on  the 
bottom  that  when  rotated  counter-clockwise  divides  the  input  by  10. 
Signal  amplification  then  is  summarized  as  follows: 


/25 

(SIG  LEV)  x(  50 
\500 


m ■ ■ ■ ■ i -■  ■ pm 

™ ® " ' ■ W ^ r5"  ' ' ■ ’ 
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(4)  When  you  are  certain  that  you  are  not  saturating  an 
amplifier,  put  TAPE  RECORDER  MONITOR  on  RECORD,  channel 
#1.  Set  POR  DIAL  to  +45  and  zero  the  meter  using  balance 
control  on  EPR  consol . 

, 

(5)  Make  a trial  run  through  entire  spectra  observing  signal 
on  scope  and  TAPE  RECORDER  MONITOR  to  make  sure  you  are 
recording  what  you  want  to  record. 

(6)  Monitor  the  control  level  on  channel  #4  to  see  if  input 

level  gives  95%  for  high  and  0%  for  low  level.  This  level 

i controls  the  PDP-8  analog  to  digital  converter  by  telling 

it  when  to  take  data. 

; 

(7)  Make  sure  channel  #2  is  on  and  internally  set  for. 

Flutter  compensation.  (Front  panel  comes  open)  This 
compensates  the  other  three  channels  for  flutter. 

(8)  Turn  off  strip  recorder, 

(9)  Put  POR  DIAL  on  -50  in  control  level  high  area. 

(10)  Punch  record  on  tape  recorder. 

(11)  Punch  increase  on  consol. 

24.  After  recording  a spectrum  go  back  to  step  n.  of  23.  and  see  if 
the  symmetry  has  shifted.  It  is  quite  fragile  and  may  be  unsatisfactory 

* 

now.  Keep  working. 

25.  Record  DPPH  at  approximately  3.32  KG  using  steps  23  e.  thru  24. 

26.  After  recording  tet  and  cubic  sites  and  DPPH  rotate  +39°  on  mag 
and  set  FIELD-SET  dials  to  1740. 

27.  Find  an  appropriate  amplification  and  identify  trigonal 
resonance.  Use  above  procedure  to  get  symmetry.  Record  trigonal  and 


kJL 


DDPH  resonance. 

28.  Go  to  Appendix  C for  transfer  of  analog  data, 
and  to  DTSS. 
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Figure  B4.  Hook-up  for  Recordi 


APPENDIX  C 

PROCEDURE  FOR  DATA  TRANSFER  TO  DTSS 


1.  Go  to  Room  E-6,  PDP-8  room,  on  Michel  son  lab  deck. 

2.  Set  up  tape  recorder,  oscilloscope,  and  PDP-8  as  shown  in 
Figure  Cl . 

3.  Refer  to  REDIFAX  #1  to  turn  on  machine  #1. 

4.  Refer  to  REDIFAX  #2  to  connect  to  DTSS  time  sharing. 

5.  You  must  have  these  programs  in  your  catalog: 

DIGACQ  DUB(EMPTY) 

CORGOS  PUB (EMPTY) 

UNPACK 

INTO  DUB 

MY  LIST  (EMPTY) 

M'Tbin  (EMPTY) 

DIGACO  and  CORGOS  are  used  on  the  PDP-8  so  they  must  be  assembled  in 
PDP-8  binary.  To  do  this  refer  to  REDIFAX  #5.  The  assembled  versions 
of  DIGACQ  and  CORGOS  can  be  saved.  I called  them  DIGACQB  and  CORGOB. 
Note: PDP-8  programs  use  octal  numbers. 

6.  DIGACQB  must  be  sent  to  PDP-8  memory.  Refer  to  REDIFAX  # 4. 
EXAMPLE 

a.  Type:  OLD  DIGACQB/NPA/PUN/PAR\0 

b.  Wait  for  a bell  and  carriage  return.  Signifying  that  trans- 
fer is  complete. 

0 s Zero. 

r Back  Slash  a Capital  L 

c.  Type:  CTRL  SHIFT  0 (Letter  0).  You  are  now  in  "ODT."  Refer 


to  REDIFAX  #4. 


I , 

i\ 
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d.  Type:  200/  Should  get:  7604  This  is  the  instruction  line 

4200  in  DIGACQB.  If  not  go  to  a.  and  start  again. 

7.  Setting  switch  register  for  Date  Sampling  Rate. 

a.  One  number  from  the  M555g  column  must  be  placed  on  the  switch 
register. 


Run  Time  M55510  M555S 


Minutes  Sec  (Base  10)  (Octal) 


1 

60 

146.43 

222 

2.5 

150 

366.21 

556 

5 

300 

732.4 

1334 

10 

600 

1464.84 

2670 

Formula  for  M555 


10 


(Run  time  (sec)  x 10^)/4096.  10^/sec  is  the  fre- 


quency of  clock  in  PDP-8  program.  4096  data  words  of  storage.  M555 


10 


must  be  converted  to  octal  M555g.  The  procedure  for  this  is  as 
fol lows: 


a.  (M555.jq/8L)  = X where  L is  largest  power  of  8 for  which 


M55510  > 3L. 


b. 

c . 

d. 


Record  the  Integer  part  of  X. 

And  multiply  the  decimal  part  of  X by  8. 

Go  to  b.  until  you  have  enough  places  of  the  octal  number. 


LX AMPLE:  146.4810  * 222, 

(1) 

(146.48]0/3Z)  = 2.28875 

Record 

2 

(2) 

(2.28875-2) (8)  = 2.31 

Record 

22 

(3) 

( .31 ) (8)  = 2.48 

Record 

222 

(4) 

( .48) (8)  = 3.34 

Record 

222.3 

(5) 

Round  to  nearest  lower  value  = 

222. 

1 


8.  Run  PDP-8  Program, 

a.  Type:  200G. 


C3 


l : 

n 


b.  Switch  register  lights  should  read  0207 

(000  010  000  111) 

c.  On  tape  recorder  put  "MONITOR"  on  "REPRODUCE  4." 

d.  Push  play  (forward)  button. 

e.  When  meter  level  goes  to  high  level  raise  switch  register  0 to 

upper  position.  (NOTE:  Scope  trace  from  Channel  #2  should  be  high 

also.)  See  257  blinking. 

f.  When  level  falls  (indicating  that  PDP-8  should  take  data)  Sw. 
Reg.  lishts  should  read  241.  You  can  monitor  data  on  scope  ch  #1. 

g.  Wait  until  377  (blinking)  comes  on,  indicating  all  data  taken, 
level  is  high  and  program  is  in  display  mode.  (000  Oil  111  111) 

h.  Detach  cable  from  scope  ch  #1  and  connect  D/A  channel  #1  to 
scope  channel  #1  and  see  data  display. 

i.  Push  Sw.  Reg.  0 down  to  see  unfiltered  data. 

j.  If  data  is  not  filtered  satisfactorily: 

(1 ) Halt  program. 

(2)  Set  7000  on  Sw.  Reg.  (Ill  000  000  000). 

(3)  Punch  ADDR  LOAD,  EXT  ADDR  LOAD,  CLEAR,  CONT.  Type:  150/ 

Should  get:  0017. 

(4)  Type:  NEW  OCTAL  NUMBER.  (See  REDIFAX  #3) 

(5)  Type:  320G 

(6)  Observe  with  Sw.  Reg.  0 = NP.  Go  to  (1)  if  not  satisfied. 

k.  If  satisfied  halt  program  and  put  4000  on  Sw.  Reg. 

l.  Punch  ADDR  LOAD,  EXT  ADDR  LOAD,  CLEAR,  CONT. 
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9.  Sending  filtered  data  to  DTSS  (See  REDIFAX  *11). 

a.  Type:  /OLD  C0RG0b/NPA/PUN/PAR\0 

b.  Wait  for  bell  and  carriage  return  from  computer. 

c.  Type:  /UNS. RESULT. /DIR 

d.  Receive:  SPEAK! 

e.  Type:  CTRL  SHIFT  0 (Letter  0) 

f.  Receive:  CARRIAGE  RETURN 

g.  Type:  2700/ 

h.  Receive:  0000 

i.  Type:  0003  for  filtered  data  or  0002  for  non-filtered  data 

j.  Type:  LINE  FEED 

k.  Receive:  2701/0000 

l.  Type:  0000  STARTING  LOCATION  OF  DATA 

m.  Type:  LINE  FEED 

n.  Receive:  2702/0000 

o.  Type:  7777  END  LOCATION  OF  DATA 

p.  Type:  CARRIAGE  RETURN 

q.  Type:  2600G 

r.  Wait  until  you  get  "READY"  (Could  be  as  long  as  5 minutes) 

s.  Type:  STR  RAN  12  (Carriage  Return) 

t.  Receive:  READY  NEW  FILE  IS  .RESULT. 

u.  Type:  /SAV/OLD  UNPACK/RUN/OLD  INTO  DUB/RUN  See  example 

printout  (Figure  C2) 

10.  You  are  now  ready  to: 


a.  Put  DUB  on  DTSS  system  mag  tape  for  storage. 

b.  Start  number  crunch  to  get  absorption  curve. 
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Ic.  Start  over  because  computer  went  down. 

11.  Output  of  DUB  is  an  APL  file,  in  the  form  DATA  = N 3.R0. 

12.  See  "BACKGROUND"  manual  for  instructions  for  magnetic  tape;  you 
must  check  out  a tape  from  Ward  Hall. 

I 
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APPENDIX  D 
DATA  CRUNCH 

1.  I assume  that  APL  data  file  "DATA  = N8.R0"  is  in  BASIC  FILE 

i 

CALLED  DUB. 

2.  Run  XAPL  (After  each  system  command  or  completion  of  program  A 
"?"  will  be  given. 

( 

* 3.  )Read  DUB. 

4.  .RODATA 

Receive:  N 3 

5.  )Read  APLLIB*** :RAWC0R 
RAWCOR  DATA 

Receive:  NEW  NAME  CORDAT 

6.  Check  first  and  last  numbers  in  0CRDAT=0 

Ex:  1 ^ CORDAT  (char  ret) 

* 0 

"1  | CORDAT  (char  ret) 

0 

If  not  look  at  data  and  go  to  5. 

7.  At  this  point  you  can  )REA  APLLIB*** :CUT  PLOT  and  look  at  data  or 

continue  with  3.  NOTE:  Result  of  CUTPLOT  is  CPOUT.  This  must  be 

written  into  a basic  file  by  ) WR1  FILE  NAME;  CPOUT.  It  may  be 

advisable  to  )SCR  FILE  MARE  first, 
c * 

8.  You  probably  need  to  write  CORDATA  into  a file  and  )EXIT.  If 

your  work  space  in  APL  is  small,  you  may  be  near  the  limit  if  you  have  4000 
PDP-8  data  numbers  in  DUB.  Then  "RUN  XAPL"  and  read  CORDATA  back 
into  work  space. 

9.  Correct  phase  distortion  in  first  derivative  data. 

a.  CUTPLOT  must  be  used  to  plot  200  of  N data  points  using 


procedure  7.  Each  distorted  first  derivative  curve  must  be  corrected 
individually.  The  scale  of  "X"  axis  plot  will  be  0 to  200.  There  are 
N/200  real  data  points  associated  with  each  point  on  the  plot.  Pick 
out  the  limits  of  your  distorted  first  derivative  curve  and  input 
these  as  limits  of  "DISTORT"  program.  Refer  to  Appendix  E for  the 
program  description.  LIMITS  = (L1)(N)/200  (L,)(N)/200 

where  L-j , are  X coordinates  from  plot. 

b.  )REA  DISTORT 

c.  DISTORT  CORDATA 

d.  Receive:  CORD  See  Step  10 

e.  CORDATA  = CORD 

f.  Repeat  c.  thru  e.  until  all  individual  distortions  are  gone. 

10.  To  make  sure  distort  worked,  see  if  .CE/CORD  = .AB.FL/CORD. 

11.  Perform  indefinite  integration  of  first  derivative  curve  to  get 
absorption  curve. 

)READ  APLLIB*** : INTTEG 
INTTEG  CORD 
Receive:  I CURV 

12.  You  may  need  to  write  ICURV  into  a file  and  exit,  run  XAPL,  read 
I CURV  back  into  work  space. 

13.  Correct  for  integration  ramp  in  absorption  curve. 

)REA  APLLIB*** :ABSC0R 

ABSCOR  ICURV 
Receive:  CORDATA 

14.  Read  7.  In  order  to  fit  a particular  absorption  curve,  it  will 
be  necessary  to  plot  the  entire  absorption  spectra  and  to  select  a 


suitable  segment  of  the  absorption  spectrum.  Follow  the  procedure 
for  determining  segment  limits  as  outlined  in  step  9.  Follow  step  7. 
to  plot  CORDATA. 

15.  Scratch  file  containing  ICURV.  Write  CORDATA  into  file  that  used 
to  contain  ICURV. 

)REA  FIT 
FIT  CORDATA 

This  will  return  the  Gaussian  parameters,  amplitude,  exponent  co- 
efficient, and  standard  deviation.  Data  is  now  in  CURFIT. 

16.  CURFIT  is  a true  Gaussian  distribution.  The  interval  of  X 
components  is  the  integer  1 so  no  special  double  precision  program 
for  integration  need  be  run, 

17.  )REA  APLLIB*** : INTSIM 
INTSIM  CURFIT 
OUTPUT  = AREA 

13.  This  completes  crunch  except  for  amplification  factors  used  in 
original  data. 


nr  ca  T » 

= True  Area . 

Amplification  Factor 


19.  Continue  with  next  data  batch  at  step  1.  or  write  paper  for 
publication  and  receive  Ph.D. 


El 

APPENDIX  E 
PROGRAMS 

DIGACQ:  A PDP-8  program  to  convert  the  analog  ESR  signal  into  digital 

numbers. 

CORGOS:  Sends  data  from  PDP-8  to  DTSS.  See  Redifax  #11. 

UNPACK:  Changes  coded  data  into  a random  access  file. 

INTODUB:  Changes  a random  access  file  into  an  APL  file. 


RAWCOR:  This  program  eliminates  the  PDP-8  modulo  409G  wrap-around, 

removes  DC  bias,  and  corrects  for  initial  potato.  Result 
is  called  CORDAT. 

CUTPLOT:  Is  used  to  cut  down  the  number  of  points  to  be  plotted. 

In  this  way  a file  of  4096  data  points  may  be  reduced  to 
200  plotting  points.  See  Appendix  D for  use  of  plotted 
graph.  Result  is  called  CPOUT. 

DISTORT:  This  program  corrects  for  phase  distortion  produced  by 

inbalance  in  phase  tuning  on  klystron.  It  assumes  a 
Gaussian  error  in  the  leakage  current  that  detects  the 
magnetic  resonance.  Result  is  called  CORD. 

INTTEG:  This  program  performs  indefinite  integration  on  the  first 

derivative  data.  The  range  of  integration  is  an  input  to 
the  program.  The  result  is  called  ICURV.  The  output  is 
the  absorption  curve  with  an  integration  ramp. 

ABSCOR:  This  program  corrects  the  integration  ramp  that  results  from 

the  indefinite  integration.  The  result  is  called  CORDATA. 


E2 


FIT: 


INTSIH: 


APLTOBAS ; 


ESR  PLOT: 


This  program  fits  a particular  absorption  curve  to  a 
Gaussian  curve,  The  boundaries  of  the  curve  are  inputs. 
Amplitude  A and  exponential  coefficient  B are  outputs-,  ST 
is  standard  deviation.  The  mid-point  of  the  Gaussian 
curve  is  found  by  the  program.  The  result  is  called  CURFIT. 
This  program  performs  the  definite  integration  on  the 
fitted  curve  to  arrive  at  a number  equal  to  the  area  under 
the  fitted  curve. 

This  program  is  necessary  for  converting  APL  data  files 
into  standard  graphic  format.  Use  CUTPLOT  while  in  APL; 
exit  from  APL;  run  APLTOBAS  (a  BASIC  program). 

This  program  outputs  the  g-splitting  or  magnetic 
resonance  point  spectra  for  CaF^iEr^*  in  table  form  or 
in  a form  suitable  for  plotting,  (a  BASIC  program).  It 
requires  "Split"  as  an  output  file,  (from  -30  to  390°) 


mm 
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DIGACQ 


I 


100  * 106 
1 10  LOOK,  0 
120  DCA  TEMPI 
130  TAD  K300I 
140  ADSC 
150  ADCV 
160  JMS  I SAMPI 
170  TAD  M2000 
180  SMA  CLA 
190  JMP  LAST1 
200  CLA 
210  TAD  K3000 
220  ADSC 
230  ADCV 
240  CLA 
250  TAD  TEMPI 
260  JMP  I LOOK 

270  LAST1 , CLA 

271  TAD  ONE 

272  TAD  STIND 

273  DCA  STIND 
280  LS.CLA 
290  TAD  I K101 
300  DCA  I STIND 
310  ISZ  STIND 
320  JMP  LS 

330  JMP  I DISPI 
340  K300I,  2001 
350  M2000,  -2000 
360  K3000,  2000 
370  M555,  -222 
380  K5400,  5400 
390  STIND,  0 
400  MATAL,  0 
410  SIG.4000 
420  L. 17 
430  K100,  100 
440  JTAL,  0 
450  TEMP,  0 
460  0 

470  TEMPI,  0 
480  K101 , 0101 
490  SAMPI,  SAMP 

5 00  DISPI, 320 

501  ONE,  -1 
510  *200 
520  LAS 
530  CIA 

540  DCA  M555 
550  LAS 


L 
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DIGACO  (continued) 

560  SMA  CLA 
570  JMP  .-2 
580  TAD  K300I 
590  A DSC 
600  ADC V 
810  JMS  SAMP 
620  TAD  M2000 
630  SPA  CLA 
640  JMP  .-3 
650  JMS  SAMP 
« 660  TAD  M200 0 

670  SMA  CLA 
680  JMP  .-3 
690  CLA 
700  OCA  STIND 
I 710  TAD  K3000 

720  A DSC 
730  ADCV 
740  STA 
750  CLZE 
760  CLA 
770  TAD  M555 
780  CLAB 
790  CLA 
800  TAD  KS400 
810  CLOE 
820  CLA 

830  6221  / CDF2 
840  LOOP  1 , CLSK 
850  JMP  .-1 
860  CLSA 
870  JMS  SAMP 
880  JMS  LOOK 
890  DCA  I STIND 
900  ISZ  STIND 
910  JMP  LOOP  I 
920  JMP  LAST) 

930  SAMP,  0 
940  ADSF 
950  JMP  .-1 
960  ADRB 
970  ADCV 
980  .IMP  I SAMP 
990  *270 
1000  DISP,  LAS 
1010  SMA 
1020  JMP  .+3 
I0?0  6231 
1040  JMP  .+2 
1050  6221 
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DIGACQ  (continued) 

1060  CLA 
1070  TAD  K3000 
1080  6552 
1090  CLA 
1100  6552 
1110  DCA  STIND 
1120  L00P2,  TAD  I STIND 
1130  6551 
1140  CLA 
1150  ISZ  STIND 
1160  JMP  L00P2 
1170  ISZ  STIND 
1180  JMP  ,-1 
1190  JMP  DISP 
1200  *320 
1210  SWAB 
1220  TAD  K100 
1230  DCA  MATAL 
1240  UPTHAR,  TAD  MATAL 
1250  DCA  STIND 
1260  TAD  L 
1270  CIA 
1280  DCA  JTAL 
1290  CAM 
1300  DST 
1310  TEMP 
1320  LOOP3,  6221 
1330  TAD  I STIND 
1340  6201 
1350  TAD  SIG 
1360  MQL 
1370  DAD 
1380  TEMP 
1390  DST 
1400  TEMP 
1410  CLA 
1420  ISZ  STIND 
1430  ISZ  JTAL 
1440  JMP  LOOP 3 
1450  CAM 
1460  DAD 
1470  TEMP 
1480  DVI 
1490  L 
1500  ACL 
1510  TAD  SIG 
1520  6231 
1530  DCA  I MATAL 
1540  6201 
1550  ISZ  MATAL 


05/14/76 


CORGOS 

100  *2600 

110  CLA 

111  TAD  K7600 

112  AND  POINT 
in  DCA  POINT 

114  TAD  K7600 

115  AND  TALLY 

116  TAD  K177 

117  DCA  TALLY 
120  TAD  FELD 
130  CLL  RTL 
140  RAL 

150  TAD  KCDF 
160  DCA  .+1 
170  CDF 
180  CLA 

190  6224  / RIF 
200  DCA  FELD 
210  LOOP.  CLA 
220  TAD  I POINT 
230  7012  / RTR 
240  7012  / RTR 
250  AND  K377 
260  JMS  SEND 
270  CLA 

280  TAD  I POINT 
290  AND  K 17 
300  CLL  RTL 
310  RTL 
320  DCA  TEMP 
330  ISZ  POINT 
340  TAD  I POINT 
350  AND  K7400 
360  CLL  RTL 
370  RTL 
380  RAL 
390  TAD  TEMP 
400  JMS  SEND 
410  CLA 

420  TAD  I POINT 

430  AND  K377 

440  JMS  SEND 

450  ISZ  POINT 

455  NOP 

460  CLA 

470  TAD  TALLY 

480  CMA 

490  TAD  POINT 

500  S2A 

510  JMP  LOOP 


05/14/76 


C0RG03  (continued) 


320  TAD  FELD 
530  TAD  KCDF 
540  DCA  ,+1 
550  CDF 

553  JMP  I K4000 

554  SEND.  0 

555  5311/  TSF3 

556  JMP  .-1 

557  6316  / TLS3 

558  CLA 

560  JMP  I SEND 
562  *2700 
564  FELD,  0 
566  POINT,  0 
568  TALLY,  0 
870  TEMP,  0 
880  M8,  -10 
890  K377.  377 
910  K7400,  7400 
020  K 17.  17 
930  K4000 , 4000 
960  KCDF,  CDF 
963  K7600,  7600 
966  K177,  177 
970  $ 
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UNPACK 

90’ THIS  PROGRAM  UNPACKS  THE  RANDOM  ACCESS  FILE  THAT  CONTAINS 

91 ’THE  DATA  FROM  THE  PDP-8  AND  RECONSTRUCTS  IT  INTO  MEANINGFUL  DATA 

92' SEE  MASS  TRANSFER  OF  DATA  IN  REDIFAX 

100  LET  1=0 

110  DIM  A(4096),S( 13) 

120  FILE  #1:  ".RESULT;" 

125  FILE  #2:  "PUB" 

130  FOR  J = 1 TO  L0F(#1) 

140  READ  #1  : S$ 

150  CHANGE  S$  TO  S 
160  FOR  K=  0 TO  3 
170  LET  1=1+1 

180  LET  A(I)=l6«S(3*K+1)+INT(S(3«K+2)/l6) 

190  LET  1=1+1 

200  LET  A(I)=256*MOD(S(3«K+2),l6)+S(3*K+3) 

210  NEXT  K 
220  NEXT  J 
230  SCRATCH  #2 
240  FOR  J=  I TO  I 
250  WRITE  #2:  A(J) 

260  NEXT  J 
270  END 
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INTODUB 

70'  THIS  PROGRAM  PUTS  4a6  ROWS  OF  8 NUMBER  DATA  INTO  FILE 
71'  DUB.  SEE  TEXT  FOR  DETAILS 

72'  TO  GET  RID  OF  AVERAGING  DEADSPACE  WE  LOSE  THE  FIRST  8 ROWS 
IV  AND  THE  LAST  8 ROWS  OF  PDP-8  DATA 
100  FILE  #1:  "PUB" 

110  FILE  #2:  "DUB" 

120  SCRATCH  A 2 

121  FOR  I = 1 TO  64 

122  READ  #1:0 

123  NEXT  I 

125  PRINT  #2:  "DATA  =";INT(LOF(#1)/8)-l6;"8.RO" 

130  FOR  N = 1 TO  LOF ( # 1 ) /8— 1 6 
140  READ  /M:X 

150  READ  #1:Y 

151  READ  #1:Z 

152  READ  #1:A 

153  READ  If  1 :B 

154  READ  # 1 :C 

155  READ  A1:D 

156  READ  #1 :E 

160  PRINT  #2:X;Y:Z:A:B;C;D;E 
170  NEXT  N 
175  PRINT  If 2 : 

180  PRINT  "DATA  IN  DUB" 

100  END 


RAWCOR 


RAWCOR  D 

.CM  ELIMINATES  PDP-8  FIELD  WRAP-AROUND,  REMOVES  DC  BIAS, 
.CM  AND  CORRECTS  FOR  INITIAL  POTATO. 


CORDAT= ,D 
D=0 

CORDAT=CORDAT+4095#-2048 . LECORDAT 
CORDAT=CORDAT- rCORDAT 
X=  1"C0RDAT 
Y="  rCORDAT 

CORDAT=CORDAT+  ( " 1 + . 10 . ROCORD AT )#(X-Y)%("1+. ROCORDAT ) 
'"CORDAT”  IS  CORRECTED  RAW  DATA  MATRIX.' 


'DIMENSIONS  OF  CORDAT  ARE' ; .ROCORDAT 


$ CUTPLOT ; BIG ; N ; NP ; NN 

[ 1 J .CM  NIL.ADIC  OPR  WITCH  SELECTS  AN  INPUTTED  NUMBER  OF 

[2]  .CM  EQUALLY  SPACED  VALUES  FROM  INPUTTED  DATA.  USED 

[3]  .CM  TO  CUT  DOWN  THE  NUMBER  OF  POINTS  TO  BE  PLOTTED. 

[4]  .CM  OUTPUT  IS  AN  N BY  2 MATRIX,  THE  FIRST  COLUMN 
15]  .CM  BEING  THE  INTEGERS  FROM  0 TO  N-1. 

[6]  •' 

1 7]  ’INPUT  DATA  VECTOR.' 

18]  BIC.r.fl 

[9]  ” 

[10]  Nr, ROBIG 

[11]  ’INPUT  NUMBER  OF  POINTS  DESIRED.’ 

[ 12]  NP=" 1+? 

[13]  NN=.FL(0i5+N*NP) 

[14]  IF  4088<NP//NN  THEN  NN=NN-1 

[15]  CPOUT= ( NP , NN ) , ROBIG 

[16]  CPOUT=CPOUT[ ; 1] 1"BIG 
[ 17]  NP=  l>NP 

[ 18]  CP0UT=.TR(2,NP).R0((" 1+. IONP) .CPOUT) 

I 19]  " 

[20]  ’PLOT  MATRIX  IS  CALLED  ’’CPOUT".’ 

[21]  NP;’  VALUES  OUT  OF’;. ROBIG;’  ARE  USED.’ 

[22]  " 

U 


DISTORT 


$ DISTORT  DAT 

[1]  ’THIS  PROGRAM  CORRECTS  FOR  PHASE  DISTORTION’ 

[2]  'INPUT  DER  CURVE  LIMITS' 

[3]  L1=0 

[4]  DAT= , DAT 

[5]  CDAT=DAT[" 1+L1[ 1 ]+i IO( 1+ — /LI ) ] 

[6]  DO=;CE/CDAT 

[7]  AO=.FL/CDAT 

[8]  XZ=" 1+L1[ 1]+CDAT.IODO 

[9]  XZN="1+L1[1]+CDAT.IOAO 

[10]  X2G=  0 i 5# ( XZ-XZN ) 

[11]  C=.CE((XZ+XZN)#O.S) 

[12]  B1=0.5*(X2G*2) 

[13]  J= ( 0 i 5#-D0+A0 ) #* ( B 1 # ( X2G ) *2) 

[ 14]  C0RD=J#*(-B1fl((-C)+;I0(.R0DAT))*2) 

[ 15]  CORD=CORD+DAT 

[16]  'DATA  IN  CORD' 


El  4 


"DISTORT'1  Inputs  the  entire  resonance  spectra  and  then  asks  for  the 
"X"  limits  between  which  you  want  to  apply  the  leakage  current  dis- 
tortion correction.  The  first  derivative  distortion  is  assumed  to  be 
caused  by  improper  phase  balance  between  phase  tuning  on  the  klystron 
and  the  energy  absorption  detection  circuits  (from  crystal  detector 
to  the  leakage  current  meter  on  the  ESR  console).  During  a run 
through  a magnetic  absorption  range  the  leakage  current  must  not 
drift  from  250yA  (an  arbitrary  constant  for  convenience).  However, 
during  a resonance,  the  leakage  current  must  increase  then  decrease 
as  the  microwave  detection  bridge  in  the  klystron  becomes  unbalanced 
indicating  that  some  energy  has  been  absorbed  as  a resonance  point 
was  passed.  Since  I am  assuming  a Gaussian  absorption  curve,  I also 
assume  that  the  leakage  current  distortion  is  an  Gaussian  curve. 

The  effect  of  this  program  is  to  simulate  the  Gaussien  leakage  current 
distortion  and  subtract  this  from  the  first  derivative  data.  See 
Figure  El  and  E2.  The  following  picture  shows  where  the 

parameters  for  di stort.  come  from  this  picture.  Note  that  |Ag|  < [Dg|. 
See  Appendix  B for  the  procedure  that  would  have  corrected  this. 


X2G  = position  for  non-distorted  distribution  first  derivative  peaks 
General  form:  y = A EXP[-B(x-C)2] 

y'  = -2BA(x-C)EXP[-B(x-C)2] 
y"  - 4B2A(x-C)2EXP[-B(x-C)2]+(-2BA)EXP[-B(x-C)2] 

At  x-C  = X2G 

y"(X2G)  --  0 = 4B2A(X^)2EXP[-B(X2G)2]  + (-2BA)EXP[-B(>(^)‘!] 
4B2A(XlG)2EXP[-B(X5G)2]  = 2BA  EX P [ - B ( X2G ) 2] 

Factor  out  2BA  EXP[-B(X2G)2]  from  both  sides 


The  amplitude  F of  the  leakage  current  Gaussian  distortion  curve  is 


proportional  to 


|An| -|D 

(|A0|-|D0I>/2 


In  this  case,  jAQ|  < |DQi,  the  correction  amplitude  would  be  negative 
which  is  correct  to  bring  first  derivative  back  to  symmetry  i.e., 
|Ag|  = i D0 I * 

It  is  not  obvious  but  the  amplitude  of  the  correction  must  be 
normalized  to  the  fitting  parameter  B.  At  X-C  = X2G 

F EXP[-B(X2G)2]  = ( I Aq | - | D Q | )/2 

F = EXP[B(X2G)2]  • ( I A0 | - | Dq l )/2 

9 

Correction  Function  = F EXP[-B(x-C)  ] 
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INTTEG 

$ INTTEG ;D;BND;S;H;N; NEXT 

[1]  .CM  NILADIC  OPER  TO  DO  AN  "INDEFINITE"  INTEGRATION  OF  A 

[2]  .CM  VECTOR  OF  Y- VALUES  FOR  EQUALLY  SPACED  X-VALUES. 

[3]  .CM  FOR  MAX  NUMBER  OF  DATA  POINTS  (SAY  4024),  MAY  USE 

[4]  ;CM  OVER  50  SECS  RUN  TIME.  AFTER  INTEGRATION  BEST  TO 

[5]  ;CM  WRITE  WHAT  YOU  NEED  INTO  FILES,  EXIT,  AND  THFN 

[6]  .CM  'RUN  APL'  AGAIN,  STARTING  FRESH  ON  RUN  TIME; 

[7]  ’’ 

[8]  ‘INPUT  VECTOR  OF  VALUES  TO  BE  INTEGRATED,' 

[9]  D=,P 

[10]  " 

[11]  'INPUT  RANGE  OF  INTEGRATION  AS  A TWO-COMPONENT  VECTOR.' 

[12]  ’THEN  WAIT  PATIENTLY.’ 

[13]  BND=0 

[14]  DD=D["1+BND[ 1]+,I01i — /BNDl 

[15]  Sr.RODD 

[16]  ICURVrS.ROO 

[17]  Hs0.5#(—  /BND)%("1+3) 

[18]  ICURV[ 1 ]=H#DD[ 1 ] 

[19]  ICURV[  2]=H#+/2”'DD 

[20]  N=2 

[21]  NEXT:N=N+1 

f 22]  ICURV[ N]  = ICURV[ N« 1 ]+H#(DD[ N- 1 ]+DD[N] ) 

[23]  ,GO(N<S)/NEXT 

[24]  ” 

[25]  'LIMITS  OF  INTEGRATION:  FROM’ ;BND[ 1 ] ; ’ TO' ;BND[2] ; • . ' 

[26]  'X-INCREMENT  IS' ;2#H 

[27]  ’OUTPUT  INTEGRATED  Y(X)  IS  CALLED  "ICURV".' 

[28]  ” 

U 


1.0 


□ 8 
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4BSC0R 


$ ABSCOR  OD 

[1]  .(>1  MONADIC  OPR  WHICH  CORRECTS  DATA  FDR  BASELINE  SLOPE, 

[?]  OD=.OD 

[3]  COPDATA=OCH.("U,IO.ROOD)#((rOD)-("rOD))t("U;ROOD) 

[4]  CORDATA=~CORDATA 

[5]  'NEW  VECTOR  CALLED  "CORDATA".' 

[6]  •• 

$K 


f 


•i 

' 

. 

i 


1 


- — 


El  9 


-1- 


FIT 

$ FIT  DAT 

[1]  .CM  THIS  PROGRAM  FITS  AN  GSR  ABSORPTION  CURVE  TO  A GAUSSIAN 

[2]  .CM  CURVE.  THE  LIMITS  BOUNDING  THE  CURVE  ARE  INPUTS. 

[33  .CM  THE  AMPLITUDE  AND  THE  EXPONENTIAL  COEFFICIENT  ARE 

[4]  .CM  OUTPUTS.  THE  STANDARD  DEVIATION  OF  THE  DATA  FROM  THE 
[53  .CM  GAUSSIAN  IS  GIVEN. 

[6]  'INPUT  PEAK  LIMITS  LI  L2‘ 

[73  L1=# 

[83  CDL1=DAT["UL1[1]+,I04~  /LI] 

[93  YMAX=.CE/CDL1 

[10]  XMAXr " 1+L 1 [ 1 3+CDL 1 ■ IOYMAX 

[11]  FD1r(CDL1>0,5#.CE/CDL1)/CDL1 

[12]  G1rDAT["1+L1[1]+iIO(XMAX-L1[1])] 

[133  FD2r(G1<0;5#.CE/CDL1) 

[14]  NPTSG1=FD2.I00 

[15]  STARTPTr ( " 1+L 1 [ 1 3+NPTSG 1 ) 

[16]  NriROFDI 

[17]  XrSTARTPT+"1+;ION 

[18]  Sr+/(X-XMAX)«2 

[19]  Q=+/.LGFD1 

[20]  Rr+/( iLGFDI )#( (X-XMAX)*2) 

[21]  Pr+/(X-XMAX)»4 

[22]  Mr 2 2iR0(N,(-S)tS,(-P)) 

[233  V=Q,R 

[24]  DATFITrM  SIMEQ  V 

[25]  X="1+L1[1]+.I0.R0CDL1 

[ 26]  DATFIT[ 1 ]r*DATFIT[ 1 ] 

[27]  CURFITrDATFIT[ l]#»-DATFIT[2]#(X-XMAX)»2 

[28]  STr(((+/(CURFIT-CDL1)»2)%("1+;R0CDL1))*0.5) 

[29]  'AMPLITUDE  =• ;DATFIT[ 1] ; • B =• ;DATFIT[2] 

[30]  'STANDARD  DEVIATION  r*;ST 

[31]  'DATA  IN  CURFIT' 

$K 

$ SOLNrCF  SIMEO  C;JDEX;DEL;M 

[1]  .CM  THIS  DYADIC  OPERATOR  OBTAINS  THE  SOLUTION  (VECTOR  VC) 

[2]  .CM  CF  A SET  OF  LINEAR  SIMULTANEOUS  EQUATIONS  — 

[33  .CM  THE  DETERMINANT  OF  ITS  COEFFICIENTS  (CF)  IS  THE 

[4]  iCM  PREOPERAND  AND  THE  VECTOR  OF  CONSTANT  TERMS  (C) 

[5]  .CM  IS  THE  POSTOPERAND. 

[6]  VCr C ; JDEXrO ; DELr DET  CF; 

[7]  NEXTJ:MrCF 

[8]  M[;JDEXrJDEX+1]rC 

[9]  VC[ JDEX]r(DET  M)IDEL 

[10]  .GO( JDEX<, ROC) /NEXT J 

[11]  SOLNrVC 
$K 

$ ZrDET  W;K 

[ 1 ] W[  K , 1 ; ]rW[  1 , KrK.  IO.CE/Kr . ABMM[ ; 1 ] ; ] 

[2]  Zr(1.E"9<.ABMM[1;1])#MM[1;1]#»1»K.NE1 

[33  .GOOfilOiOR/ ( 1 .EQ.ROMM) ,O.EOZ 
[4]  ZrZ#DET  1 1.DAMM-MM[;1].SO.#MM[1;]*MM[1;1] 

*K 


"FIT."  Inputs  the  entire  magnetic  absorption  resonance  spectrum. 

You  can  then  specify  the  X coordinate  limits  which  contain  the  absorp 
tion  curve  to  be  fitted  to  a Gaussian  curve.  It  then  finds  all  the 
points  in  the  set  greater  than  half  the  height  of  the  peak.  The 
amplitude  (A)  and  exponential  coefficient  (B)  are  found.  The  mid- 
point of  the  curve  is  taken  as  the  average  X coordinate  distance 
between  the  negative  and  positive  first  derivative  peaks.  Therefore 
M is  known  from  data.  The  following  is  the  derivation  for  A and  B. 

Gaussian  Curve  Fit: 

y = A EXP[-B(x-m)2] 

M = Midpoint  of  the  Gaussian  curve;  found  from  data 
In  y = In  A+(-B)(x-m)2 
F = C+(-B)(x-m)2 
Di  = ith  Data  Point 

Least  Squares  Sum: 

z(F-ln  Di)2  = S(C,B) 


i 
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- — = z[C-B(x-m)2-Z]  = 0 
2 3C 

\ .. 

— ■ l[C-B(x-m)2-Z][x-m]2  = 0 
2 3B 

zC  - Bz(x-m)2  - zZ  = 0 
Cz(x-m)2  - Bz(x-m)4  - zZ(x-m)2  = 0 

Let  z(x-m)2  = S 

Z (Z) (x-m)2  = R 
z(x-m)4  = P 
ZC  = Czl  = CN 

CN  - BS  - Z = 0 

CS  - BP  - R = 0 

|n  - si  r ci  fzi  r ci  fN  - sr1  r z] 

Is  - pj  IbJ  LrJ  [b]  Is  - pj  Lr J 

A = EXP(C)  = EXP(ln  A) 
Therefore  A and  B are  found. 
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* INTSIM  DAT 

[.5]  .CM  THIS  PROGRAM  PERFORMS  DEFINITE  INTEGRATION 
[ • 6 J .CM  ON  INPUT  DATA  FILE  USING  SIMPSON'S 
t .73  .CM  RULE.  THE  PROGRAM  FINDS  AN  EVEN  NUMBER 
[.81  .CM  OF  SEGEMENTS  BY  DROPPING  THE  LAST  DATA  POINT 

[1]  S=6 
[1.1]  DATArDAT 

[2]  X= , 102 

[5]  .G0( ( .RODAT)  ,NE  2#.FL( ( .RODAT) )*?)/NEXT 

[6]  DATA=" 1;DADAT 

[7]  NEXT : N= .FL( "0 . 5+ • RODATA ) 

[9]  I=(  (" l+. RODATA )t( 3#N) )#+/( 1 , ( (N-1) ,R04  2),1)#DATA 

[ 10]  'THE  COMPUTED  VALUE  OF  THE  INTEGRAL  IS:' 

[ 11]  S SIGNIF  I 

$K 

$ SIG=S  SIGNIF  B;A;C;M;N 

[1]  N=IV2#B<0 

[2]  A=N#(B+B,EQO) 

[ 3]  C= .FL10.LG, A 

[4]  M=( 10.LG, A)-C 

[5]  SIG=(.ROA).RO(.FL(0.5+10»M+S-1))#10*1+C-S 

[6]  SIG=N# ( SIG-B . EOO ) 

$K 
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APLTOBAS 

100'  THIS  PROGRAM  CHANGES  ANY  SET  OF  fN  2;R0]  APL  DATA  MATRICES 
110'  IN  A GIVEN  FILE  INTO  STANDARD  GRAPHIC  FORMAT; 

120'  WRITTEN  BY  DANNY  MILLS  *76  (MAR  1976). 

130  DIM  V(100) 

140  PRINT  "APL  DATA  FILE  NAME"; 

1*0  LINPUT  X$ 

160  FILE  #1:  X$ 

170  FILE  #2: 

180  LINPUT  #1:  Y$ 

190  LINPUT  #1:Y$ 

200  IF  Y$  <>  " " THEN  230 
210  PRINT  #2:  "1E37.1E37" 

220  GO  TO  410 

230  CHANGE  Y$  TO  V 

240  FOR  K1  = 1 TO  V(0) 

250  IF  V(K1 ) <>  34  THEN  270 
260  LET  V(K1 ) = 45 
270  NEXT  K1 

280  FOR  K1  = 1 TO  V(0) 

285  IF  V(K1)=45  THEN  310 
290  IF  V(K1)  = 32  THEN  310 
300  GO  TO  320 
310  NEXT  Kf 

320  FOR  LI  = K1  TO  V(0) 

330  IF  V(L1)  = 32  THEN  345 

335  IF  V(L1 )<>  45  THEN  360 

336  CHANGE  V TO  V$ 

337  LET  Yi  = SEG$(V$,1,L1-1) 

338  LET  Z*  = SEG$(V$,L1,LEN(V$)) 

339  LET  V$=  Y$4CHR$(44)4CHR$(32)AZ$ 

340  CHANGE  V$  TO  V 

341  GO  TO  370 
345  LET  V(L1)=44 
350  GO  TO  370 
360  NEXT  LI 

370  CHANGE  V TO  Y$ 

380  PRINT  #2:  Y$ 

390  IF  END  #1  THEN  430 

400  GO  TO  190 

410  IF  END  #1  THEN  430 

420  GO  TO  180 

430  RESET  #2 

440  SCRATCH  #1 

450  FOR  I = 1 TO  7000 

460  LINPUT  #2:  M 

470  PRINT  #1:  R$ 

480  IF  END  #2  THEN  500 
490  NEXT  I 

500  PRINT  "CHANGED  TO  STANDARD  GRAPHIC  FORMAT" 

520  END 
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ESRPLOT 

100  LET  H =.7 147788*9. 348 

105  LET  0 = 1 

110  FILE  #1:  "SPLIT'' 

116  SCRATCH  #1 
120  LET  SI  = 6.25 
130  LET  S2  = 7.78 
140  LET  T1  = 8.55 
150  LET  T2  = 3.29 
160  LET  P = 3.14159/180 

165  PRINT  " TABLE  OR  GRAPH"; 

166  INPUT  X$ 

167  IF  X$="GRAPH"  THEN  170 

168  IF  X$="TABLE"  THEN  500 

169  GO  TO  165 

170  PRINT  "B  OR  G" ; 

171  INPUT  Y$ 

172  IF  Y$="B"  THEN  179 

173  LET  H = 1 

174  LET  Os-1 

179  FOR  T=-30  TO  390  STEP  5 

180  LET  S3  = SQR( (S2*C0S( P*T) )~2+(S1*SIN( P*T) )*2) 

190  PRINT  #1:  T;",";H/S3~Q 

200  NEXT  T 

210  PRINT  01:  "1E37,1E37" 

220  FOR  T = -30  TO  390  STEP  5 

230  LET  S4  = S0R(S2“2«.5*(SIN(P*T))“2  ♦ Sr2»(1-  i5*(SIN(P*T))"2)) 
240  PRINT  #1:  T;",";H/S4“0 
250  NEXT  T 

260  PRINT  01:  "1E37.1E37" 

270  FOR  T = -30  TO  3go  STEP  5 
280  LET  B = (T-54.75)*P 

290  LET  S5  = SOR( (T2»C0S(B) )“2  + (T1»SIN(R) )~2) 

300  PRINT  01:  T;",";H/S5“Q 
310  NEXT  T 

320  PRINT  01:  "1E37,1E37" 

330  FOR  T = -30  TO  390  STEP  5 
340  LET  C r (T+54.75)»P 

350  LET  S6  = SQR( (T2»C0S(C) )"2  ♦ (T1»SIN(C) )~2) 

360  PRINT  01:  T;",";H/S6“Q 
370  NEXT  T 

380  PRINT  01:  "1E37.1E37" 

390  FOR  T s -30  TO  390  STEP  5 
400  LET  C = (T+54.75)»P 

410  LFT  S7  = SOR(  1/3»(T2*COS(P«T))“?  VlT2»(  1-1/3«(C0S(  P*T) )"2) ) 
420  PRINT  01:  T:",";H/S5“Q 
430  NEXT  T 

440  PRINT  01:  "1E37.1E37" 

450  FOR  T = -30  TO  390  STEP  5 

460  LET  S8  = 6.78 

470  PRINT  01:  T;"." ;H/S8*0 


ESRPLOT  (continued) 


480  NEXT  T 

490  PRINT  #1:  "1E37.1E37" 

499  GO  TO  700 

500  PRINT  "B  OR  G"; 

501  INPUT  Y$ 

502  IF  Y$="B"  THEN  510 

503  LET  H = 1 

504  LET  Q=-1 
5 10' CONTINUE 

511  PRINT  #1:  " T TET1  TET2  TRIG1  TRIG2  TRIG3  CUBIC" 

512  FOR  T=  -30  TO  390  STEP  5 

521  LET  S3  = SQR ( ( S2*C0S( P*T ) ) * 2+ ( S 1 *SIN ( P*T ) ) * 2 ) 

531  LET  S4  = SQR(S2“2*»5*(SIN(P*T))~2  + Sr2»(1-  i5*(SIN(P*T))“2)) 

541  LET  B = (T-54.75)*P 

542  LET  S5  = SQR((T2«C0S(B))“2  + (T1*SIN(B))“2) 

551  LET  C = (T+54.75)»P 

552  LET  S6  = S0R((T2»C0S(C))*2  + (T1*SIN(C))“2) 

562  LET  S7  = S0R(1/3»(T2*C0S(P»T))“2  +Tr2»(  1-1/3*(COS(P»T))“2) ) 

581  LET  S8  = 6;78 

582  LET  S3=INT((H/S3)"0»10000+;5)/10000 

584  LET  S4=INT((H/S4)*Q*10000+;5)/ 10000 

585  LET  S5=INT((H/S5)A0*10000+;5)/10000 

586  LET  S6=INT((H/S6)*Q»10000+;5)/10000 

587  LET  S7=INT((H/S7rQ*10000+.5)/10r  ' 

588  LET  S8=INT( (H/S8)*Q*10000+»5)/1u^ou 

599  PRINT  #1:  T;S3;S4;S5;S6;S7;S8 

600  NEXT  T 
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APPENDIX  G 

DERIVATION  OF  TRANSFER  FUNCTION 
OF  DIGITAL  MOVING  AVERAGE  FILTER 


The  digital  moving  average  filter  is  basically  a moving  window 
memory  L locations  wide.  It  sumsthe  contents  of  the  first  L memory 
locations  in  field  2,  divides  by  L,  and  stores  this  number  in  the 
Lth  memory  location  in  field  3 of  the  PDP-8  memory.  The  window  then 
moves  one  number  to  the  right  and  runs  the  contents  of  memory  loca- 
tions 2 thru  l+l , divides  by  L,  and  stores  the  result  in  the  L+lth 
memory  location.  This  continues  through  the  entire  memory.  The 
idea  is  to  make  the  window  contain  an  integral  number  of  wavelengths 
of  the  signal  to  be  filtered  out.  The  action  of  the  filter  in  the 
time  domain  can  be  written  in  terms  of  the  Kth  memory  location.  Each 
location  holds  the  value  of  some  time  varying  function  at  time  (K)(S) 
where  K is  the  number  of  the  sample  (e.g.  the  5th  sample  taken:  K=  ) 
and  S is  the  time  between  samples  taken  by  the  analog  to  digital 
converter  on  the  PDP-8  minicomputer.  (1/S  is  the  sampling  frequency.) 
For  L samples  [Kth  to  (K-L+l)-th]. 

Y(K)  = 1-[X(K)+X(K-l)+X(K-2)  + ...X(K-L+2)+X(K-L+l)] 

L 

Now  take  the  Z- transform  of  both  slides 


Y(z)  = 1-[X(z)Z0+X(z)Z'1  + ...+X(z)Z"L+1] 


= X.(jl(i+z"1+z"2+..I+Z_L+1) 

L 


The  transfer  function  = Y(z)/X(z)  = H(z) 


H(z)  = 


loz‘1+z“2+..,z"L+1) 

L 


Using  the  identity 


K=L-1 
I Z 
K=0 


-K  _ 1-Z 


-L 


1-Z 


„(z)  z 

L 1-Z  IZL  Z-l  Lr"1  Z-l 


This  is  visualized  in  the  Z domain  as  follows: 
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Let  Z = eJ<1) 


H(jiu)  = 


1 ejuL 


[e'j«a^T][e’j».r] 


jwL  jwL  jwL 

_ e TJ  e 7 -e^7  ]_ 

Jw  _ jw 

L[eJ“(L‘5)][e2  - e 2] 


juL  jwL 

i hu)  i ■ 

Jw  Ju 

J _“T, 


e - e 


Using  Euler's  identities 


eJw  = cosu  + jsinu 
e'Ju  = cosu  - jsinu 


|H(u) | 


sin(!f) 
L sin(|) 


For  u=0  j H ( 0 ) | = 

sin(O) 


Using  L' Hospital's  rule 


I H(qj)  l..=n  = Lini 


Leo 

sin  -o- 


Z.„Umki  c.°lr 


(*>■►0  L sin(^)  w-*T)  L/2  cos  ^ 


For  w^O  and  |H(u)|  = 0 sin  ^ =0  or  ^ = 


Therefore  u * f = n 


t 
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if  ,the  Nyquist  frequency. 


For  an  unstable  condition  (pole  in  denominator) 

sin^=0  or  ^=mr  therefore  go  = 2nir  , 

but  at  go  = 2nn  the  numerator  equals 

sin^ip-  = sin  mrL  = 0 

for  all  integers  n and  L.  As  stated  before 

sin(O)  , 
sin(O) 


There  are  no  unstable  frequencies  for  the  filter.  The  frequency  that 
is  filtered  out  is  given  by 


U> 


2 Iff  = GO 


f = 


n 

L 


I chose  L = 15  (17  in  octal)  so  that  any  data  periodic  over  15 
samples  is  filtered  out.  The  sampling  time  was 


10x60 

mr 


.146 


seconds  or  a frequency  of  6,3/sec. 

The  Myquist  frequency  is  then  3.4  Hz.  A typical  resonance  of  a 
large  peak  takes  30  seconds  with  the  range  set  at  1000  Guass  and  sweep 
time  set  for  10  minutes.  This  gives  a frequency  of  .033  Hz  = f and 
go  = 2uf  = ,066lT. 
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tables  of  b and  g 

ROAD  MAP  VALUES 


SPlITTIilG  VALUES 

-1- 


112 


AS  A FUNCTION  OF  "T" 
ANGLE  OF  ROTATION 


T TET1  TET2  TRIG  I TRIG2  TRIG3  CUBIC 
-30  7.4271  6.4611  8i 5195  4.6626  7;585  6.78 

-25  7.5296  6,4015  8.4339  5.1146  7; 4867  6.78 
-20  7.6169  6;  3496  8.2942  5i573  7.-4008  6.78 

-15  7.687  6.3073  8.102  6.024  7; 3303  6.78 
-10  7.7383  6.2758  7.8594  6i4569  7.278  6,78 

-5  7.7695  6.2565  7.5695  6;8631  7.2457  6.78 

0 7.78  6.25  7;2359  7.2359  7i2348  6.78 
5 7.7695  6.2565  6.8631  7.5695  7.2457  6;78 


10 

7.7383 

6.2758 

6.4569 

7.8594 

7.278 

6*78 

15 

7.687  ( 

>.3073  < 

i;024  8 

i 102  7.: 

3303  6; 

78 

20 

7.6169 

6.3496 

5.573  1 

3.2942  ‘ 

7 ; 4008 

6.78 

25 

7.5296 

6.4015 

5.1146 

8.4339 

7 i 4867 

6.78 

30 

7.4271 

6i 46 1 1 

4.6626 

8.5195 

7;  585 

6.78 

35 

7.3121 

6;5264 

4.235  1 

3.5499  1 

7.6924 

6i78 

40 

7.1874 

6.5952 

3; 855  i 

8.525  7 

.8051  6 

.78 

45 

7.0566 

6.6655 

3;5511 

8.4449 

7.9198 

6i?8 

50 

6.9233 

6.735  : 

3543  l 

8,3106  < 

3.0328 

6.78 

55 

6;7915 

6i 80 18 

3.2902 

8; 1235 

8.1408 

6,78 

60 

6.6655 

6i8638 

3.3683 

7.8859 

8.2409 

6.78 

65 

6.5495 

6.9194 

3; 5772 

7.6005 

8.2303 

6.78 

70 

6.4478 

6.9671 

3;8901 

7.2711 

8.4068 

6;78 

75 

6.364  ' 

r.0055  * 

1.2761  1 

6i902  8 

*4683  6 

.78 

80 

6.3016 

7.0336 

4.707  1 

6.4989  1 

3*5133 

6.78 

85 

6.263  1 

N0508  1 

j. 1604  1 

6,0683  1 

3.5408 

6;78 

90 

6.25  7. 

,0566  5, 

,6186  5 

.6186  8 

.55  6.78 

95 

6i?63  1 

f.0508  ( 

j.0683  1 

5;  1604  i 

3,5408 

6.78 

100 

6.3016 

7;0336 

6.4988 

4.7071 

8.5133 

6.78 

105 

6.364 

7.0055 

6.902 

4.2761  1 

3.4683 

6*78 

110 

6.4478 

6.9671 

7.2711 

3; 8901 

8.4068 

6,78 

115 

6.5495 

6.9194 

7.6005 

3.5772 

8.3304 

6.78 

120 

6.6655 

6.8638 

7.8858 

3.3683 

8;2409 

6.78 

125 

6.7915 

6.8018 

8.1235 

3.2902 

8,1408 

6i78 

130 

6.9233 

6;735 

8.3106 

3.3543 

8.0328 

6;78 

135 

7.0566 

6.6655 

8.4449 

3.5511 

7.9198 

6*78 

140 

7.1874 

6.5952 

8.525 

8.855  ' 

7.8052 

6.78 

145 

7.3121 

6.5264 

8.5499 

4;  235 

7.6924 

6.78 

150 

7.4271 

6,4611 

8.5195 

4.6626 

7.585 

6*78 

155 

7.5296 

6.4015 

8.4339 

5.1146 

7i4867 

6.78 

160 

7.6169 

6 i 3496 

8.2942 

5.573 

7.4008 

6.78 

165 

7.687 

6.3073 

8.102 

6.024  7 

.3303  6k78 

170 

7.7383 

6.2758 

7.8594 

6.4569 

7;  278 

6.78 

175 

7.7695 

6.2565 

7.5695 

6.3631 

7.2457 

6.78 

180 

7.78  1 

3.25  7.2359  7. 

2359  7i 

2348  6. 

78 

185 

7.7695 

6; 2565 

6.8632 

7.5695 

7.2457 

6i78 

190 

7.7383 

6; 2758 

6.4569 

7.8594 

7.278 

6*78 

1 95 

7.687 

6. 3073 

6.024 

8.1019 

7.3303 

6.78 

200 

7.6169 

6.3496 

5.573 

8.2942 

7 ; 4008 

6.78 

205 

7.5296 

6,4015 

5.1146 

8,4339 

7.4867 

6.78 

210 

7.4271 

6.4611 

4.6626 

8.5194 

7.585 

6;78 

■■MM 
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H3 


SPLIT  (continued) 


215  7.3121  6.5264  4i235  8;5499  7.6924  6i78 

220  7.1874  6.5952  3.855  8.525  7i805 1 6.78 

225  7 i 0566  6.6655  3.5511  8.4449  7*9198  6*78 

230  6.9233  6*735  3;3543  8.3106  8.0328  6*78 

235  6.7915  6.8018  3;2902  8.1235  8.1408  6.78 

240  6.6655  6.8638  3i 3683  7.885Q  8*2409  6i78 

245  6.5495  6,9194  3*5772  7*6005  8.3303  6*78 

250  6.4478  6.9671  3.8901  7*2711  8.4068  6.78 

255  6i364  7*0055  4.2761  6.902  8,4683  6i78 

260  6.3016  7.0336  4.707  6.4989  8.5133  6*78 

265  6*263  7.0508  5*1604  6;0683  8.5408  6.78 

270  6.25  7.0566  5*6186  5i6l86  8.55  6*78 
275  6,263  7.0508  6.0682  5.1604  8*5408  6*78 

280  6.3016  7*0336  6*4988  4.7071  8.5133  6i78 

285  6; 364  7.0055  6*902  4.2761  8.4683  6.78 

290  6.4477  6.9671  7*2711  3.8901  8*4068  6.78 

295  6.5495  6.9194  7.6005  3*577 2 8.3304  6'.78 

300  6i6655  6.8638  7.8858  3.^683  8;2409  6;78 

305  6.7915  6;8018  8.1235  3.2902  8i1408  6i78 

310  6.9233  6; 735 1 8.3106  3.3543  8.0328  6;78 

315  7.0566  6.6655  8;4449  3.5511  7.9198  6i78 

320  7il874  6.5952  8.525  3;855  7.8052  6.78 

325  7.3121  6.5264  8.5499  4i235  7.6924  6.78 

330  7.4271  6i 46 1 1 8.5195  4;6626  7i5851  6.78 

335  7.5296  6.4015  8;4339  5.1146  7j4867  6.78 

340  7.6169  6.3496  8.2942  5i573  7i4008  6.78 

345  7i687  6.3073  8,102  6.024  7;3303  6i78 
350  7.7383  6;2758  7.8594  6.4569  7.278  6;78 

355  7.7695  6.2565  7i5695  6.8631  7.2457  6;78 

360  7;78  6.25  7.2359  7;2358  7.2348  6,78 
365  7.7695  6;2565  6.8632  7.5695  7;2457  6.78 

370  7.7383  6.2758  6i4569  7.8594  7i278  6;78 
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EE424  - REDIFAX  #1 


TURNING  MACHINE  ON  6 OFF : STARTING  6 STOPPING  A PROGRAM 


'I he  Key  (on  a chain)  switch  controls  the  power. 
Turn  to  mid  position  only  for  normal  use. 

Turn  to  left  for  off. 


To  start  a program  running: 

Turn  machine  on. 

Set  starting  address  on  switen  regisrer. 
Press  I.OAi)  AL)DR . 

Press  CLEAR,  CONTINUE. 


Usual  starting  address  are: 


0200  - Working  program 
4000  - 8SIIARE 
7000  - odt 


See  REDIFAX  #463 


in  stop  a program  simply  press  HALT.  A halted  program  may  be  restarted 
(rom  beginning  as  above  or  allowed  to  continue  from  where  it  was  by 
pressing  CONTINUE. 


lo  cause  a program  to  go  one  instruction  at  a time  keep  HALT  switch 
down  and  press  CONTINUE  repeatedly. 


- Be  sure  to  restore  the  system  before  shutting  down.  (See  REDIFAX  6 ). 


lor  more  details  on  starting  programs  (e.g.  starting  a program  in  a different 
Field)  see  'Small  Computer  Handbook'. 


0 
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EE424  - REDIFAX  #2 


CONNECTING  TO  TIME  SHARE  SYSTEM 


- Machine  I 


Turn  on  '.Jefferson  Box'  and  punch  in  A,  D,  F,  G,  H,  J,  § K.  (Pressing 
M clears  any  former  key  in) . 

Turn  on  Beehive  Terminal  (push  button)  and  press  Half  Duplex  Button. 
The  Beehive  is  optional  but  provides  a useful  monitor. 

Dial  up  TSS  using  left  hand  modem  and  usual  procedure  (Dial  8 only). 

If  8SHARE  is  not  already  running,  start  it  (See  REDIFAX  4 ). 

Communication  with  TSS  should  be  evidenced  by  sign-on  message. 


- Machine  II 


Turn  on  the  machine  and  start  8SHARE  (REDIFAX  1&2). 

Dial  up  TSS  using  right  hand  modem  and  usual  procedure  (Dial  8 only). 
Communications  should  be  evidenced  by  sign-on  message. 

Both  Machines  - Pressing  TALK-CLEAR  on  related  modem  will  disconnect. 

Connection  remains  intact  even  if  machines  are  not  in  8SHARE  and 
communication  may  be  resumed  by  simply  starting  8SHARE.  Of  course, 
if  no  activity  occurs  over  a long  period  TSS  will  disconnect.  Restart 
as  above. 


EE424  - RED I FAX  #3 


0DT  HIGHLIGHTS 


- To  start  0DT  see  RELIIFAX  1.  Start  location  is  7000.  Normally  resident 
in  fieldO. 

- To  0pen  Register  ABCD  : type  A BCG/  , where  ABCD  is  the  octal  address. 

0DT  responds  with  the  octal  cont  i.ts. 

- To  Modify  an  0pen  Register  : type  the  desired  new  octal  value  followed 
by  a carriage  return  or  line  feed. 

- To  Close  an  0pen  Register  without  Modification  : type  carriage  return 
or  line  feed. 


- To  0pen  Next  Register  in  Sequence  : type  line-feed.  This  will  also  close 
the  previously  opened  register  (with  or  without  modification).  0DT 
responds  with  octal  address/octal  contents  of  next  register. 

- To  start  a program  at  octal  location  WXYZ,  type  WXYZG.  0DT  gives  a 
line-feed/carriage  return  and  passes  control  to  stated  location. 

- To  place  a breakpoint  at  octal  location  WXYZ  type  WXYZB-  0DT  responds 
with  a line-feed/carriage  return  and  retains  control.  When  you  start 
your  program  running  the  instruction  at  WXYZ  will  be  replaced  by  a trap 
which  will  return  control  to  0DT  when  program  control  reads  WXYZ 
(instruction  at  WXYZ  has  not  been  executed). 

- To  continue  running  from  a breakpoint,  type  0C.  0DT  responds  with  line- 
feed/carriage return,  executes  the  original  instruction  at  the  break- 
point , and  transfers  control  back  to  the  program. 


- io  remove  a breakpoint,  type  OB. 


To  transfer  control  from  01)T  to  8SHARE, 


type  C0NTROI.-SHIFT-0  (all  at  once). 


- To  transfer  control  from  8SHARE  to  0DT, 


type  CONTROL-SHIFT - L (all  at  once). 


i 


■ 


For  complete  details  on  0l)T  see  Chapter  5,  pp  68-83,  of  Introduction 
To  Programming 
and 

EE  Dept  document  '8SHARE1  which  details  local  modifications  to  0DT. 


"0" 

"0" 


represents  zero. 

represents  the  alphabetic  character  0. 


NOTH : 


EE424  - REDIFAX  #4 


8SI1ARE  HIGHLIGHTS 


To  start  8SHARE  and  connect  to  TSS  see  REDIL'AX  1 IJ  2.  Starting  location 
is  4000.  (Will  operate  in  fieldO  only). 

To  dialogue  with  TSS  .‘Start  8SHARL  and  connect  to  TSS.  Intera'ction  from 
keyboard  is  almost  the  same  as  - 7 you  were  using  any  other  ASR33/TTY  in 
the  yard. 

Exceptions : 

(a)  Bells  from  TSS  don't  ring  and  in  general  TSS  control  characters 
are  ignored. 

(b)  Lower  case  text  from  TSS  prints  out  garbage. 

(c)  Back  slash  (upper  case  L)  and  back  arrow  (Upper  case  0)  are 

reserved:  be  very  careful  to  avoid  these  unless  you  intend  a 

binary  data  transfer  - see  below. 

(d)  Only  about  one  page  may  be  listed  without  restarting,  because 
of  speed  differentials  and  limited  buffering. 

(e)  The  only  way  to  stop  transmission  from  TSS  is  by  typing  S : 
the  other  break  mechanisms  available  at  most  yard  terminals 
don't  work  with  8SHARE. 

(f)  C0NTROL-SHIFT-0  passes  control  to  0DT . 

To  transfer  a binary  program  from  TSS,  follow  this  example:  (Assuming 

dialogue  with  TSS  is  already  in  effect). 


Parity  - On . 
Off  I 

{ ^Punch 
/0I.1)  MYBIN/NPA/PUN/PAR4 

t t_i 

inary  file  j 

or  transfer  f 


J/ 

R\0 


Target  Data  Field 


\ 


Back  Slash 


Bi 

fo 


Slashes 


8SHARE  will  respond  with  a CR/LF  after  the  Target  Data  Field  is  typed 
and  the  keyboard  will  go  dead  until  the  transfer  of  the  binary  file  is 
complete.  (The  bells  that  ring  just  before  the  keyboard  reactivates 
may  be  i gnored)  . 

To  transfer  a binary  file  to  TSS,  see  8S11ARE  document. 


lor  full  details  on  8SHARF.  see  the  EE  Dept  Document  '8SHARE'. 
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8 MAP  HIGHLIGHTS 


- To  assemble  a program  using  8.MAP  you  must  have  (on  your  TSS  account): 

(a)  A source  file  of  the  mnenomic  form  instructions  you  wish  to 
have  assembled. 

(b)  A saved  file  which  will  receive  the  assembly  listing. 

(c)  A saved  file  which  wil,  receive  the  binary  program. 

Basically  8MAP  takes  the  source  file,  checks  its  consistency,  and 
produces : 

(a)  a statement  of  errors  found. 

(b)  a resultant  binary  program  in  a form  that  may  be  loaded  into 
the  PDP-8 . 

(c)  an  assembly  listing  of  the  binary  program  which  shows  the  source 
code,  and  the  way  it  has  been  assembled  (i.e.  location  and  value 
of  each  instruction). 

- To  use  8MAP  type,  RUN  8MAP . 

It  will  respond  with: 

SOURCE  ; 1. 1 ST  I NG  ; B I NARY  ; ERROR;  SYMTAB  - - 

At  this  point  you  are  expected  to  supply  at  least  three  pieces  of 
information,  i.e.,  the  file  names  for  your  source  file,  listing  file, 

L binary  file,  in  that  order.  Usually  your  current  file  will  be  the 
source  file  - indicate  this  by  a ; and  then  type  the  name  of  your  listing 
and  binary  file  separated  by  a ; and  followed  by  a CR.  (Ignore  the  ERROR 
and  SYMTAB  files). 

8MAP  then  takes  over,  prints  any  error  messages  and  builds  the  listing 
H binary  files,  and  returns  control  to  you  with  the  same  file  current. 
Previous  contents  of  files  specified  for  listing  and  binary  are  lost. 

Full  details  on  8MAP  may  be  obtained  by  typing  /0LD  PDP8LIB*** .-D0CB/LIS . 
(Be  prepared  for  a 6 or  7 page  print-out). 

Illustrative  Example: 

To  create  a 'do  nothing'  prog '-am. 


f 


-Opening  Required  Files,  if  not  already  available. 


/NtW  JdNQTc/bA V/NEW  hYLlb/lSAV/NEW  MYblN/SAV/BRl 
JSPLD  DUNOTE't 

*J00  *200  l 

110  JMH  . ( Building  Source  Program 

120  $ J 


/REH/kUN  bllAH  4- 


_Calling  8MAP 


Current  File  implied  for  SOURCE 


SOJhCtjLlbTlNb;  ul  NARY  ; EHhOR  ;SYMTAo--:  MYLISs  MYbl  N 

NORMAL  TE^hfNATION  ' 


- Underscore  signify  material  typed  by  TSS 


a 

a 


EE424  - RED  I FAX  #6 


RESTORING  THE  SYSTEM  SOFTWARE: 

(See  also  Crash  Recovery  - REDIFAX  #7) 

MACHINE  I:  The  operating  system  for  Machine  I is  on  a tape  cassette 

labelled  MASTER  SYSTEMTAPE.  Locate  and  mount  this  cassette 
on  either  drive,  if  not  already  mounted.  Basic  steps  are: 

- Press  Rewind  (unlabeiied  black  button  beside  drive) 

- Set  4540  on  'switch  register  (SR)  * 

- Press  ADDR  L0AD 

- Set  SR  to  0000 

- Press  EXTD  ADDR  L0AD 

- Set  SR  to  0100  if  cassette  on  drive  B,  otherwise  leave 
SR  at  0000. 

- Press  CLEAR/CONTINUE. 

Tape  should  read  in  and  stop  with  4600  showing  on  top  set  of 
lights  - if  not  see  REDIFAX  #7.  Any  user  programs  previously 
resident  will  be  lost. 

MACHINE  II:  The  operating  system  is  on  a full-folded  paper  tape  labelled 

MASTER  SYSTEM  TAPE-MACHINE  II.  Locate  and  mount  this  tape 
in  the  high  speed  reader  with  leader  over  the  read  sensors. 
(Leader  has  only  the  outside  hole  punched). 

The  basic  steps  are: 

- Set  7777  on  SR. 

- Press  ADDR  LOAD. 

- Set  SR  to  0010. 

- Press  EXTD  ADDR  LOAD. 

- Press  CLEAR/CONTINUE. 

lull  details  on  the  use  of  the  Cassette  Handler  are  contained  in  the 
LE  Department  document  "DUMPLOAD". 


EE424  - RED I FAX  #7 


CRASH  RECOVERY  PROCEDURES: 

Various  crash  status' may  occur  depending  on  the  machine  and  the 
extent  of  the  problem.  The  basic  idea  is  that  if  the  loader  necessary 
for  recovery  to  the  next  higher  level  is  intact,  then  there  is  no  need 
to  go  further  down  the  recovery  procedure. 

- MACHINE  I: 

Status  DI:  Simple  loss  of  system  but  with  DUMPLOAD  intact.  Follow 

procedure  given  in  REDIFAX  #6. 

Status  Cl:  Loss  of  system  including  DUMPLOAD,  but  with  BIN  LOADER 

intact . Locate  the  binary  tape  DUMPLOAD  and  read  in  with 
procedure  given  in  REDIFAX  #8.  Machine  will  then  be  in 
state  DI . 


Status  BI:  Loss  of  system  including  DUMPLOAD  and  BIN  LOADER,  but  with 

RIM  LOADER  intact.  Locate  the  tape  BINARY  LOADER  and  read  in 
with  the  procedure  given  in  REDIFAX  #9.  Machine  will  then  be 
in  status  Cl. 

Status  AI:  Total  crash  : Key  in  RIM  LOADER  as  described  in  REDIFAX  #9. 

Machine  will  then  be  in  status  BI. 

- MACHINE  II: 


Status  Cl I : Loss  of  system  but  BIN  LOADER  intact  in  Field  I.  Follow 

procedure  given  in  REDIFAX  #6. 

Status  BI 1 : System  and  Field  I BIN  LOADER  lost.  Locate  and  read  in 

the  tape  BINARY  LOADER  as  per  REDIFAX  #9.  Machine  will  then 
be  in  state  Cl  I . 


Status  All:  Total  crash.  Key  in  RIM  LOADER  as  described  in  REDIFAX 

"9.  Machine  will  then  be  in  status  B 1 1 . 


RE424  - RED I FAX  #8 


USING  THE  BINARY  LOADER: 


Lo-Speed 


(This  description  applies  to  all  cases  except  the  read-in  of 
Machine  II's  system  tape  for  which  REDIFAX  #6  applies). 


- Place  the  binary  tape  in  the  low  speed  reader  with  leader 
(Only  outside  hole  punched)  on  the  read  holes. 

- Press  START  on  reader. 

- Set  SR  to  0000  and  press  EXTD  ADDR  LOAD. 

- Set  SR  to  7777  and  press  ADDR  LOAD. 

- Press  CLEAR  S CONTINUE. 

Paper  Tape  should  read  in,  and  halt  in  the  trailer. 


Hi -Speed  - (Machine  II  only) 

- Press  binary  tape  in  hi-speed  reader  with  leader  over  the 
read  holes. 

- Set  SR  to  0010  and  press  EXTD  ADDR  LOAD. 

- Set  SR  to  7777  and  press  ADDR  LOAD. 

- Set  SR  to  3777. 

- Press  CLEAR  5 CONTINUE. 


Paper  Tape  should  read  in  and  halt  in  the  trailer. 


Note:  For  full  details  on  use  of  binary  loader  see  'Introduction  to 

Programming ' . 


EE424  - RF.DIFAX  #9 


ENTERING  AND  USING  THE  RIM  LOADER: 


4 


The  RIM  LOADER  is  the  'last  line  of  defense'  in  the  sense  that  it  is  a 
very  simple  program  which  we  can  key  in  to  read  in  a special  format  tape 
for  a more  powerful  program  (BIN  LOADER). 

Reading  in  the  BIN  LOADER  tape  wi*  . RIM: 

- Locate  the  tape  BINARY  1.0ADER  and  mount  on  the  low  speed  reader  with 
leader  (only  outside  hole  punched)  over  the  read  holes. 

- Set  SR  to  0000  for  Machine  I or  0011  for  Machine  II. 

- Press  EXTD  ADDR  LOAD. 

- Set  SR  to  7756  and  press  LOAD  ADDR. 

* Set  the  low  speed  reader  to  START. 

- Press  CLEAR  f,  CONTINUE 

Tape  should  read  in  but  does  not  come  to  a stop  on  its  own. 

Entering  RIM  1.0ADER: 


Should  the  misfortune  of  a total  system  crash  befall  you: 

- Set  SR  to  0000  for  Machine  I or  0011  for  Machine  II  and  press  EXTD 
ADDR  LOAD . 

- Set  SR  to  775b  and  press  ADDR  LOAD. 

- Locate  the  plaque  on  the  machine  face  labelled  RIM  LOADER  and  enter  the 
octal  code  shown  in  the  LOW  column  as  follows: 

(a)  Set  SR  to  given  octal  code. 

(b)  Lift  DEPOSIT  switch. 

(c)  Repeat  at  (c)  until  all  code  is  entered. 

(d)  Check  code  by  setting  SR  to  7756  and  pressing  LOAD  ADDR. 

(e)  Set  rotary  switch  to  Ml). 

(f)  Press  EXAMINE;  lower  lights  will  display  code  in  location  one 
before  location  displayed  on  upper  lights. 

(g)  Repeat  (f)  until  all  instructions  are  checked. 

(h)  If  a discrepancy  is  found  it  is  recommended  that  you  start  afresh. 

(i)  Read  in  the  BINARY  LOADER  as  detailed  above. 
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17  December  1975 
REDIFAX  #11 


From:  Major  D.  A.  Wright 

To:  All  users  of  EE  Dept  Minicomputers 

Subj:  Mass  data  transfers  to  time  share 

1.  A new  system  program  is  now  available  which  facilitates  large  data 
transfers  to  TSS.  This  new  program  "C0RG0"  has  two  advantages  over 
transfers  using  8 SHARE.  First,  it  moves  data  33%  faster;  second,  the 
resultant  data  file  can  be  easily  converted  to  a form  amenable  to 
processing  under  BASIC,  by  means  of  a subsidary  program  "UNPACK." 

(These  programs  are  available  under  the  mini-lab's  account  number 
E02100  and  are  both  saved  public.) 

2.  Prior  to  moving  data  C0RG0  must  be  brouaht  down  to  the  PDP8  by 
8 SHARE,  (i .e./0LD*E021 OO:C0RG0/NPA/PUN/PAR~Q) . When  ready  to  move 
the  data  a target  file  named  "PUB"  must  be  opened  and  readied, 

(i.e./OLD  PUB/SCR/DIR).  Time  share  will  respond  with  (3)  SPEAK!  and 
wait  for  the  data  transfer  to  begin. 

3.  The  parameters  of  the  data  transfer  must  be  conveyed  to  C0RG0  as 

follows:  Go  to  0DT,  (i.e.,  press  CONTROL-SHIFT  0)  and  enter  the  field, 

start  location  and  end  location  for  the  data  transfer  in  locations 
2700-2702  respectively.  The  actual  data  transfer  is  then  initiated 

by  passing  control  to  2600:  (i.e.,  2600G) , and  proceeds  at  the  rate 

of  20  PDP-8  words/second. 

4.  When  the  data  transfer  is  complete  C0RG0  will  pass  control  to 
8 SHARE  (recognizable  by  a carriage  return)  and  after  a 10  second 
time-out,  time  share  will  respond  "READY."  Don't  touch  the  keyboard 
during  time-out,  as  anything  typed  will  be  appended  to  the  transfer. 

5.  The  resultant  file  (.RESULT.)  should  be  converted  to  Random  Access 
Numeric  format  by  the  following  steps: 

STR  RAN  12 

/SAV/0LD  UNPACK/RUN. 

Upon  completion  the  file  PUB  will  contain  one  entry  for  each  word  of 
core  tranferred  and  will  be  unsigned  (i.e.,  4000  thru  7777  octal  will 
show  up  as  2048  thru  4095).  If  signed  numbers  are  desired,  the  sign 
can  be  easily  restored  in  any  program  which  processes  the  file. 

6.  If  you  want  to  preserve  the  data  long-term,  be  sure  to  get  it  into 
another  file,  else  the  next  usage  of  PUB  will  destroy  it. 

7.  Note  that  C0RG0  only  passes  full  pages  of  core,  regardless  of  the 
start  and  end  locations.  That  is,  the  transfer  starts  at  the  first 
location  of  the  page  containing  the  start  address  and  ends  on  the  last 
location  of  the  page  containing  the  end  address. 
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